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ABSTRACT 


We  present  deterministic  and  probabilistic  models  for  the  analysis  of  strategic  strikes 
against  transportation  networks.  The  deterministic  models  use  integer  programming  to  solve 
problems  on  single  and  multicommodity  networks.  The  aims  of  a  network  interdictor  are 
(a)  to  minimize  the  maximum  network  flow  with  a  fixed  amount  of  interdiction  resources, 
or  (b)  to  minimize  the  total  effort  and  mission  turnaround  time  if  given  sufficient  resources 
to  stop  the  flow  completely.  In  the  case  of  a  multicommodity  network,  the  interdictor  also 
aims  to  utilize  minimum  resources  to  achieve  a  disconnecting  set  which  severs  the  paths 
connecting  all  sources  to  their  respective  sinks.  In  the  probabilistic  model,  arc  capacity  is 
not  a  factor  and  the  objective  of  a  single  interdictor  is  to  minimize  the  probability  of 
infiltration  by  a  single  evader  through  a  network  while  the  objective  of  the  evader  is  just  the 
opposite. 
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I.  INTRODUCTION 


A  strategic  interdiction  campaign  against  the  logistics  supply  system  of  an 
enemy  may  take  many  forms;  it  may  consist  of  attacks  on  the  key  installations  such 
as  military  industrial  bases,  supply  depots,  or  attacks  on  the  vital  links  between  these 
installations  and  the  front-line  troops  in  the  battlefield.  The  key  installations  are 
normally  well  protected  through  hardening  and  the  deployment  of  air  defenses  and 
troops.  Compared  to  the  installations,  the  land  transportation  network,  as  it  traverses 
through  a  large  expense  of  territories,  is  harder  to  protect  against  concerted  and 
determined  efforts  by  the  air,  naval  and  land  forces. 

A.  OBJECTIVE 

The  objective  of  this  thesis  is  to  develop  models  for  the  optimum  allocation  of 
multiple  types  of  strike  resources  to  interdict  a  land  transportation  network  in  order 
to  degrade  the  performance  of  an  enemy’s  logistics  supply  system.  In  particular,  we 
will  be  formulating  mathematical  programs  for  the  strikes  against  single  commodity 
and  multicommodity  networks. 

B.  OVERVIEW 

Strategic  strike  is  a  complex  task  for  the  Armed  Services  to  undertake.  It 
involves  the  setting  of  campaign  goals,  development  of  attack  plans  and  coordinated 
execution  by  the  participating  strike  forces. 
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1.  T^'pes  of  Strike  assets 

The  resources  available  to  the  operational  planners  at  the  command 
headquarters  for  the  interdiction  campaign  consist  of: 

a.  Air  power  such  as  a  fleet  of  fighter  aircraft  and  attack  helicopters, 

b.  Naval  gunfire  support, 

c.  Special  forces  trained  for  deep  strike,  sabotage  and  reconnaissance 
missions,  and, 

d.  Long  range  missiles  launched  from  land,  air  and  sea  platforms. 

In  developing  a  campaign  plan,  the  strike  planner  would  estimate  the  weapons  efforts 
required  to  interdict  the  targets  based  on  a  desired  confidence  level  of  mission 
success.  The  available  strike  assets  should  then  matched  in  an  optimum  manner  to 
the  targets  according  to  the  weapon  suitability  and  economy  of  efforts,  which  in  turn 
depend  on  a  host  of  factors  such  as  the  ease  of  target  detection  and  identification, 
enemy  defenses,  target  vulnerability,  the  potential  collateral  damages. 

2.  Campaign  Objectives 

The  central  task  of  the  strike  planner  is  the  optimum  allocation  of  the 
multi-dimensional  forces  to  interdict  the  following  types  of  networks: 

a.  Single  Commodity  Network.  To  degrade  the  performance  of  a  road 
network,  a  viable  strategy  is  one  of  counter-capacity.  The  aim  is  to 
minimize  the  throughput  of  network  by  allocating  a  fbced  amount  of 
resources  to  interdict  the  most  "profitable"  arcs.  Another  approach 
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is  to  minimize  the  demands  of  the  front-line  units,  or  the  demands 
prioritized  according  to  geographical  locations.  Yet  another  type  of 
single  commodity  network  is  one  faced  in  an  unconventional  warfare 
setting.  Here,  the  enemy  attempts  to  infiltrate  through  the  network 
by  avoiding  detection;  arc  capacity  is  not  a  factor.  The  aim  of  the 
strike  planner  is  to  minimize  the  probability  of  enemy  infiltration, 
b.  Multicommodity  Network.  The  attacker  would  want  to  reduce  as 
much  as  possible  the  weighted  flow  of  multiple  commodities  from 
certain  nodes  in  the  supply  network  to  other  nodes. 

C.  LITERATURE  REVIEW 

The  network  interdiction  problem  has  been  studied  for  quite  some  time  now. 
There  are  two  common  approaches  to  modelling  the  problem;  either  a  deterministic 
model  or  stochastic  /  game  theoretic  model. 

Mustin  and  McMasters  1967  develop  an  algorithm  for  optimally  interdicting  a 
transportation  network  using  limited  assets.  They  consider  a  model  in  which  the 
capacity  of  an  arc  (/,;)  can  be  reduced  by  one  unit  by  the  expenditure  of  units  of 
resources.  The  objective  is  to  minimize  the  maximum  flow  obtainable  between  two 
nodes  s  and  t  subject  to  the  consumption  of  no  more  than  E  units  of  resources. 
Their  methodology  requires  that  the  network  be  planar  so  that  the  topological  dual 
can  be  taken  for  the  enumeration  of  many  cut-sets.  This  is  an  undesirable 
requirement.  The  source-sink  planar  graph  requires  not  only  planarity  but  also 
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requires  that  the  source  and  the  sink  lie  on  the  outer  face  of  the  graph,  which  is 
quite  a  strong  assumption.  Furthermore,  the  algorithm  is  not  easy  to  generalize  to 
handle  multiple  resources  with  which  to  interdict  arcs,  or  multiple  resources  which 
must  be  applied  together  to  interdict  an  arc. 

Nugent  and  McMasters  1969  develop  an  embellished  version  of  the  above 
model  in  that  the  damage  function  is  now  assumed  to  be  exponential.  Associated 
with  each  arc  (/,/),  there  is  a  measure  of  vulnerability  b,j.  The  proportional  reduction 
in  the  arc  capacity  with  units  of  resources  allocated  is  {1  -  exp(-b,ye,y)}.  The  object 
is  to  maximize  the  reduction  in  the  maximum  flow.  The  algorithm  also  requires  the 
topological  dual  of  the  network  to  be  taken,  so  that  the  "shortest"  routes,  which 
represent  minimum  capacity  cut-sets  in  the  primal  network,  can  be  determined. 
Since  the  damage  function  is  now  exponential,  the  allocation  of  weapons  in  one  cut¬ 
set  is  determined  using  a  Lagrangian  method.  The  algorithm  always  produces 
non-integer  solutions.  This  model  suffers  the  same  weaknesses  as  the  previous 
model.  The  method  of  taking  a  topological  dual  and  enumerating  cut-sets  is 
inefficient  if  the  network  becomes  too  large. 

Wollmer  1970  presents  two  algorithms  for  targeting  strikes  against  a  network. 
The  enemy  is  assumed  to  have  the  policy  of  either  maximizing  flow  between  a  source 
node  5  and  a  sink  node  t  or  meeting  a  given  flow  at  minimum  cost.  For  his  first 
algorithm,  the  cost  of  traversing  arc  (/,/)  is  assumed  to  be  a  linear  function  of  flow. 
After  the  arc  has  been  struck,  there  is  a  resulting  cost  increase  per  unit  time  Ac  of 
a  minimum  cost  circulation  flow  and  repair  time  t,y.  Associated  with  each  arc,  there 
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is  a  strike  value  v,^  which  is  defined  as  the  repair  cost  kjj  plus  fyAc.  The  model  is  a 
one-strike  algorithm  in  that,  at  any  one  time,  only  one  arc  which  is  that  arc  with  the 
maximum  strike  value  will  be  interdicted.  Multiple  strikes  are  approximated  by 
repeated  applications  of  the  same  algorithm.  The  second  algorithm  treats  the  arc 
costs  as  piece-wise  linear  functions  with  one  break  point.  Wollmer’s  algorithms  do 
not  produce  optimal  results  in  a  multiple  strike  scenario. 

Capps  and  Taylor  1970  develop  a  model  similar  to  the  above  model.  The 
damage  function  is  now  assumed  to  be  exponential.  Against  an  arc  (i,j)  with  a 
vulnerability  measure  bjp  the  proportional  reduction  in  capacity  d^j  due  to  units  of 
effort  is  (l-exp{hy€y}).  The  arc  repair  time  and  repair  cost  are  assumed  in  turn 
to  be  linear  functions  of  d^j.  As  in  Wollmer’s  model,  this  model  selects  and 
interdicts  the  arc  with  the  maximum  strike  value  which  is  the  sum  of  kjj  and  r,yAc. 
The  model  allows  for  arc  repairs  and  allocates  the  resources  on  a  daily  basis.  This 
model  suffers  the  same  weaknesses  as  the  previous  model. 

Preston  and  Howard  1970  develop  a  procedure  for  determining  the  optimum 
allocation  of  aircraft  against  a  transportation  network  system  using  the  technique  of 
dynamic  programming.  An  exponential  relationship  between  arc  capacity  and 
interdiction  effort  is  assumed.  The  model  requires  a  topological  dual  of  the  network 
to  be  constructed  and  a  range  of  cut-sets  enumerated.  The  minimum  capacity  for 
each  level  of  available  resources  E  =  1,2,3,...  is  determined  by  dynamic  programming. 
The  algorithm  will  increase  the  value  of  £  until  either  the  maximum  \  .me  of  £  is 
reached  or  when  the  benefit  to  be  gained  from  assigning  the  last  aircraft  is  less  than 
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the  cost  of  assigning  the  aircraft.  The  model  entails  the  enumeration  of  many  cut¬ 
sets  and  does  so  for  a  large  number  of  E.  It  is  also  rather  difficult  to  produce  for  a 
particular  war-time  scenario  a  reasonable  estimate  on  the  cost  of  assigning  one 
aircraft  and  the  benefit  of  each  mission. 

Helmbold  1971  develops  a  model  for  the  allocation  of  limited  strike  assets 
among  the  arcs  of  a  transportation  network  by  analyzing  its  topological  dual.  The 
damage  function  is  assumed  to  be  a  step  function.  A  particular  arc  (/,/)  would  have 
its  capacity  reduced  to  a  lower  level  of  capacity  C2,y  by  an  interdiction  effort  e,y. 
The  optimum  allocation  of  strike  assets  to  the  arcs  in  a  cut-set  is  solved  by  dynamic¬ 
programming.  He  claims  that  the  algorithm  can  be  extended  to  handle  the  problem 
of  multiple  resources,  but  this  would  be  possible  only  in  a  limited  way.  Furthermore, 
the  model  is  again  restricted  to  planar  networks. 

Ghare,  Montgomery  and  Turner  1971  dc  /elop  a  branch  and  bound  algorithm 
for  the  interdiction  of  a  set  of  arcs  in  order  to  minimize  the  maximum  flow  of  a 
network  subject  to  fixed  amount  of  resources.  For  each  arc  (/,;)  with  capacity  Cy,  its 
destruction  requires  units  of  resource;  no  partial  destruction  of  an  arc  is  allowed. 
The  model  computes  the  lower  bounds  on  the  maximum  flow  once  an  arc  is 
destroyed,  which  affects  the  decision  whether  to  include  or  exclude  a  particular  arc 
from  the  branching  procedure.  This  choice  of  lower  bound  for  the  branch  and  bound 
procedure  is  not  ideal.  As  a  whole,  this  model  is  based  on  a  rather  weak 
characterization  of  an  interdicted  network,  and  would  be  difficult  to  implement,  for 
instance,  with  multiple  strike  resources. 
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Wollmer  1968  develops  a  non-deterministic  model  for  determining  where  to 
place  interdicting  forces  in  order  to  maximize  the  probability  of  preventing  an 
opposing  force  (the  evader)  from  proceeding  from  one  node  in  a  network  to  another. 
The  interdictor  determines  a  probability  ity  of  placing  a  single  interceptor  on  arc  (/,/) 
given  a  known  probability  of  of  stopping  a  single  evader  traversing  the  arc  given 
that  the  evader  is  traversing  the  arc.  The  aim  of  the  single  evader  is,  knowing  the 
interdictor’s  strategy,  to  pick  a  path  which  maximizes  the  probability  of  successful 
traversal  through  the  network.  For  a  small  network,  this  problem  can  be  analyzed 
in  terms  of  a  two  person  zero-sum  game  which  can  be  easily  solved  by  linear 
programming.  For  a  large  network,  the  game  matrix  becomes  too  large  to  handle. 
Wollmer  therefore  proposes  a  solution  technique  which  is  based  on  marginal  analysis 
or  the  steepest  ascent  approach.  Initially,  all  values  of  7r,y  are  set  to  zero  and  are 
incremented  in  such  a  way  that  the  ratio  of  the  increase  in  Tt,^  and  the  decrease  in  K 
is  minimized.  Without  any  formal  proof,  he  claims  that  the  solution  is  exact  for  a 
scenario  of  a  single  interceptor  versus  a  single  evader  and  yields  an  approximately 
optimal  answer  for  multiple  interceptors.  The  main  weakness  of  the  model  is  its 
inability  to  handle  multiple  infiltrators  traversing  the  network  simultaneously.  The 
optimal  stratcgN’  of  the  evader  is  also  not  given  in  the  model. 

Danskin  1962  develops  a  game  theory  model  of  convoy  routing.  The  analysis 
solves  a  two-person  zero-sum  game  in  which  one  player,  the  convoy  command,  directs 
two  routes;  the  first  represented  by  a  vector  x  for  the  merchant  vessels  and  the 
second  by  a  vector  y  for  escorts.  The  other  player,  the  enemy  submarine  command, 
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controls  one  vector  z.  The  payoff  P,  in  terms  of  tonnage  delivered  /  dollar  cost,  is 
defined  as  the  ratio  of  T  the  total  tonnage  delivered  on  a  set  of  routes  and  C  the 
total  route  cost  including  the  cost  of  losses.  The  convoy  commander  seeks  to 
maximize  P  which  is  a  function  of  (x,  y,  2)  whilst  the  submarine  commander  seeks 
to  minimize  P.  The  model  determines  a  critical  turnaround  time  using  a 
Lagrangian  method.  If  a  route  i  has  its  turnaround  time  a,  >  Oj,  it  is  not  used.  If  a, 
<  Uo,  it  is  used.  The  main  difficulty  of  this  model,  when  applied  to  network 
interdiction  problem,  is  that  the  determination  of  Og  and  other  threshold  values  would 
require  the  solving  of  equations  involving  the  sum  of  non-linear  functions  of  T,  and 

C,  over  all  paths  through  the  network. 

D.  MODEL  ASSUMPTIONS 

The  assumptions  made  for  the  development  of  the  various  model  throughout 
this  paper  are  as  follows: 

a.  Damage  function.  It  is  assumed  that  there  would  be  no  partial 
destruction  of  a  target.  The  capacity  of  an  arc  will  be  reduced  to  zero 
upon  interdiction,  which  is  a  reasonable  assumption  with  the  advent  of 
"smart"  weapons.  Other  types  of  damage  functions  may  approximated  by 
expanding  arcs  into  sets  of  independent  parallel  arcs. 

b.  Time  frame.  We  are  assuming  that  once  the  link  is  destroyed,  the  arc 
capacity  is  reduced  instantaneously  to  zero.  The  enemy  could  employ 
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other  support  means,  such  as  constructing  pontoon  bridges  next  to  the 
destroyed  bridges,  but  this  would  take  some  time, 
c.  Independent  Arc.  It  is  assumed  that  the  destruction  of  one  arc  has  no 
impact  on  the  performance  of  other  arcs  in  the  network, 
e.  Node  Capacities.  We  assume  that  each  node  will  have  infinite  capacity. 
If  interdiction  of  a  node,  say  a  road  intersection,  is  a  possibility,  this  can 
be  accomplished  by  splitting  the  node  i  into  two  nodes  ij  and  /,  and 
connecting  a  capacitated  arc  between  these  two  nodes.  In  the  case  of  a 
directed  network,  all  arcs  originally  entering  node  i  now  enter  node  ij  and 
all  arcs  leaving  node  i  now  leave  The  model  can  also  be  generalized 
to  undirected  networks. 
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II.  DETERMINISTIC  MODELS  -  SINGLE  COMMODITY  NETWORKS 


In  this  chapter,  integer  programs  are  formulated  for  the  interdiction  of  a  single 
commodity  network.  The  aims  of  the  attacker  are  to  reduce  to  the  maximum  flow 
of  materials,  minimize  the  resource  used  to  cut  off  all  supplies  to  the  enemy  and 
minimize  the  maximum  flow  of  supplies  to  the  units,  prioritized  according  to  their 
tactical  values.  Most  of  the  models  discussed  here  are  generalizations  of  those 
described  in  Wood  1991. 

A.  NETWORK  FLOW  BASICS 

Let  G  =  {N,A)  be  a  directed  graph  with  node  set  N  and  arc  set  >1.  An  arc  is 
an  ordered  pair  (/,/)  where  /,/  i  N.  Associated  with  each  arc  (/,/),  there  is  a  flow 
capacity  u,^.  The  lower  bound  of  arc  capacity  is  assumed  to  be  zero.  Let  f  be  the 
amount  of  flow  in  the  network  from  the  source  node  s  to  sink  node  t.  The  maximum 
flow  problem  can  be  stated  as: 
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Maximize 


/ 


(1.1) 


Subject  to 


'/ 

if  /  =5 

0 

if  ^  5  or  r 

if  /  =  r 

V  (/,;•) 

O 

Al 

V  (/J) 

where  the  sums  and  inequalities  are  taken  over  existing  arcs  in  the  network. 

Let  S  be  any  set  of  nodes  in  the  network  such  that  S  contains  node  s  but  not 
t.  letT  =  N  -  S.  Then  (S,T)  s  {(ij)  :i  e  S,j  e  T}  isa  cut-set  separating  node  s  from 
node  t. 

The  dual  of  the  maximum  flow  problem  is  as  follows: 

Minimize  T,  T,  (1-2) 

•  J 

Subject  to  o,  -  =  1 

a,  -  >  0  V  (/,/) 

>  0  V  (/,/) 

o,  unrestricted  V  i 

where  a  is  the  dual  variable  corresponding  to  the  conservation  equations  and  is  the 
dual  variable  corresponding  to  x  <  u. 
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In  the  single  commodity  flow  problem,  it  is  well  known  from  the  Maximum- 
flow/minimum-cut  theorem  by  Ford  and  Fulkerson  that,  for  a  given  maximum  flow 
/,  a  cut  (S,T)  is  minimal  if,  and  only  if  all  the  arcs  of  (S,T)  are  saturated  while  the 
arcs  in  the  form  (T,S)  are  flowless  with  respect  to/.  Hence,  we  are  guaranteed  an 
integer  optimal  solution  if  the  arc  capacities  are  integer.  This  also  arises  from  the 
fact  that  the  constraint  matrix  of  the  maximum  flow  problem  is  unimodular.  We 
shall  now  develop  a  constructive  proof  for  the  following  lemma; 

LEMMA  1  The  dual  of  the  maximum  flow  problem  has  a  solution  in  which  all  dual 
variables  are  0  or  1. 

Assume  that  G  is  connected  and  that  we  have  solved  the  maximum  flow 
problem  using  a  bounded-variable  simplex  algorithm  as  in  Bradley,  Brown  and 
Graves  1977.  However,  instead  of  starting  w-ith  a  full  artificial  basis,  assume  that  full 
rank  has  been  achieved  by  discarding  the  last  flow-balance  constraint  in  the 
maximum  flow  problem.  This  implies  that  Qf^i  =  0.  Now,  any  basis  which  is 
associated  with  the  flow  balance  constraints  from  a  network  flow  problem  can  be  put 
into  the  upper  triangular  form.  Let  B  be  the  optimal  basis  in  triangulated  form  and 
assume  that  variables  have  been  reflected  about  their  upper  bounds  as  necessary  to 
put  all  +  Is  on  the  diagonal  of  B.  Reflection  does  not  change  the  dual  solution.  Let 
a  denote  the  dual  row  vector  of  length  INI  -  1  and  let  Cg  denote  the  row  vector  of 
length  IM  -  1  containing  the  costs  associated  with  the  optimal  basis.  If  =  0  it 
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follows  that  a  =  0.  If  not,  Cg  =  (000...010...0)  where  the  1  is  in  position  /  say,  and 
a  can  be  computed  from 

OtB  =  Cg 

Now  =  1,  and  we  let  p(/)  <  /  denote  the  row  of  the  off-diagonal  element  in 
column  /  of  B.  If  there  is  no  off-diagonal  element  p(i)  =  0  and  we  define  Oq  s  0. 
The  above  equation  can  now  be  solved  by  substitution; 

ay  =  0 

■  ^p(2)  =  0  =9  a,  =  0  since  a,,  =  ay  =  0  and  p(2)  =  0  or  1 


a,.y  -  a^yy.yy  =  0  =»  =  0  sloce  a^  =...=  a,._,  =  0  and  p(l-l)  <1-2 

a,  -  a^yy,  =  1  =>  a,  =  1  since  a^  =...=  a,.y  =  0  and  p(l)  <  /  -  1 
®/  +  /  ■  ^p(l*l)  =  0  »  ay+y  =  1  if  =  1  or  »  ay  +  y  =  0  if  Op(l  +  })  =  0 


®/v*7  ■  ^p(Mi)  -  0  ^  tt/vi;  -  1  if  <^p(Mi)  -  1  or  =»  ay^^y  -  0  if  a^^/^yy^  -  0 
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Thus,  a,  =  0  or  a,  =  1  for  all  /  e  N.  Next,  for  a  specified  arc  (/,;)  note  that  (a)  > 

0,  (b)  fiij  appears  in  only  one  inequality  in  the  dual  formulation  (1.2),  and  (c)  the 
objective  coefficient  associated  with  fijj  is  positive.  It  therefore  follows  that 

1.  if  a,  =  0  and  a,  =  0,  then  =  0, 

2.  if  tt,  =  1  and  oij  =  0,  then  =  0, 

3.  if  a,  =  0  and  aj  =  1,  then  =  1,  and 

4.  if  a,  =  1  and  Oj  =  1,  then  =  0. 

Thus,  the  dual  of  the  maximum  flow  problem  has  an  optimum  solution  in  which  all 
dual  variables  are  0  or  1.  QED 
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B.  MINIMIZE  MAXIMUM  FLOW 


1.  Problem  Statement 

The  logistics  planner  of  an  enemy  wishes  to  force  through  as  much  of  a 
single  commodity  from  a  source  node  to  a  sink  node  in  a  directed  and  capacitated 
network.  Given  insufficient  resources  to  completely  cut  the  supply  network,  the 
problem  of  the  interdictor  will  be  to  minimize  the  maximum  amount  of  flow  that  the 
enemy  can  force  through  the  unbroken  arcs.  The  strike  planner  can  allocate  any  one 
type  of  resource  available  to  interdict  any  one  arc.  Simultaneous  strikes  by  multiple 
resources  on  a  single  target  are  not  necessary.  (This  model  can  be  generalized  to 
include  a  co-ordinated  attack  of  multiple  strike  resources  against  a  single  target.) 
Destruction  of  an  arc  is  always  assumed  to  be  total. 

2.  Model  1  -  Network  with  Directed  Arcs 

For  an  interdicted  network  model,  the  indices  used  are: 

i,j  :  nodes, 

/  :  the  type  of  strike  resources  available,  e.g.,  cruise  missiles,  aircraft, 

special  forces,  etc. 

The  data  of  the  model  are: 

e^j,  :  amount  resource  of  type  /  required  to  destroy  arc  (/,;) 

Ej  :  the  total  strike  assets  available  for  resource  type  /, 

:  the  capacity  of  arc  (i,j). 
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The  decision  variables  are; 

7yy  :  the  binary  variable  which  is  1  if  arc  (/,;)  is  destroyed  by  resource  type 
/,  and  is  otherwise  0, 

Xij  ;  flow  through  arc  (i,j), 

/  :  total  flow  through  the  network. 

The  Min-Max  problem  that  solves  the  interdictor’s  problem  is  as  follow’s: 

Minimize  Maximize  /  (2.1) 

Y «  r 


The  second  constraint  of  (2.1)  ensures  that  the  capacity  of  arc  (/,/)  is  reduced 
to  zero  once  7,^,  is  equal  to  one.  The  solution  of  the  Min-Max  problem  will  yield  an 
optimum  solution  corresponding  to  the  maximum  flow  given  that  the  total  resources 
E  are  insufficient  to  reduce  the  capacity  of  any  cut-set  (5,7)  in  the  network  to  zero. 
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To  solve  the  above  Min-Max  problem,  we  can  take  the  dual  of  the  minimization 
problem  and  convert  the  Min-Max  problem  into  a  non-linear  program  as  follows: 


Minimize  Minimize 
Y  e  r 


Z  E  u.  (1  - 


I  (<J) 


(2.2) 


Subject  to 


a  -  a  =  1 


a.  -  or,  +  >  0 


V  (/,/) 


>-  0 

tt.  unrestricted 


V  (/J)  ,  I 

V  I 


r  =  {  y,ji  I  'I',;/  ^  (0^1)  V  (^./)  ’  ^  .  and  E  e,.,  7,,/  <  E,\l  1  } 


(‘J) 


ijl  Ujt 


where  a  is  the  dual  variable  associated  with  the  flow  conservation  constraint  and  j3 
is  the  dual  variable  associated  with  the  second  constraint  of  (2.1). 

Now,  let 

”  (1  ”  ijl)  ^  ijl 

Substituting  the  above  into  the  second  constraint  of  the  (2.2),  we  have, 

^  ^  '^^.jl^.jl  ^  0 

f  / 
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TTie  second  constraint  of  (2.2)  is  now  non-linear.  However,  we  can  prove  thaii 


or  , 

=0  (2.3) 


From  Lemma  1,  we  know  that  for  fixed  values  of  7,^/  the  optimum  solution  of 
the  inner  minimization  problem  has  equal  to  0  or  1.  This  is  true  since,  for  fixed 
7y7,  problem  (2.1)  is  just  a  maximum  flow  problem  with  multiple  upper  bounds  on 
each  arc  and  such  a  problem  can  be  solved  by  throwing  out  all  bounds  but  the 
smallest.  (Throwing  out  a  bound  sets  some  to  0.)  Let  p  be  the  event  {  7,^  =  0}, 
and  q  be  the  event  {  =  1}.  Equation  (2.3)  is  true  if  p  or  ^  is  true.  It  is  false,  if 

and  only  if  both  p  and  q  are  false,  i.e.,  only  when  {  7,^/  =  1}  and  {  =  0}.  Looking 

back  at  the  objective  function  of  (2,2),  if  {  fi/jj  =  0}  we  may  assume  that  {  7,^;  =  0} 
because  setting  7,^/  =  1  would  not  affect  the  value  of  the  objective  function  but  it 
would  unnecessarily  use  up  e,j,  units  of  resources.  Therefore,  w-e  may  linearize  the 
second  constraint  of  (2.2)  by  replacing  the  non-linear  term  7^/^^;  with  7,^^  and  the 
constraint  becomes: 

/  ■  / 
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We  now  have  the  following  mixed  integer  program  (MIP); 


Minimize 

E  E 

/  (v) 

Subject  to 

OH  -  a  =  1 
t  s 

a.  -  tt.  +  E/j 

^  1 

a,  *  ^  0 

V  (i,/) 

{>•]) 

V  / 

o 

Al 

V  (/,;•)  .  / 

f  (0^1) 

V  (/,/)  <  I 

a,  unrestricted 

V  / 

(2.4) 


Since,  for  this  model,  we  will  select  at  the  most  one  type  of  resource  to  interdict 
one  particular  arc,  we  have  for  each  arc  (/,/).  T.jj-  <  1.  In  an  arc  (/,/)  in  a  cut-set, 
there  may  be  at  most  1  of  lij,  terms  taking  on  the  value  of  1.  Hence,  D,  h^j,  can  be 
replaced  by  a  single  term 
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The  final  formulation  becomes; 


Minimize  E  u,,h^.  (2.5) 

(w) 

Subject  to  a,  -  =  1 

®,  -  *  '".j  *  r7,j,  5  0  V  (i,/) 

E  ^  V  / 

o.j) 

h.j>0  M{  ') 

7y/e(0,l)  V  (/,/),  / 

a.  unrestricted  V  / 


We  can  interpret  the  above  model  as  being  based  on  a  modified  dual  of  the 
maximum  flow  linear  programming  formulation.  Essentially,  an  (S,T)  cut  has  been 
identified  with  all  a,  =  1  for  /  c  7  and  a,  =  0  for  i  e  S.  The  value  of  E,7,y,  =  1  if  (/,/) 
is  a  forward  arc  across  the  cut  and  that  arc  is  to  be  broken;  hjj  is  1  if  {i,j)  is  a  forward 
arc  across  the  cut  but  that  arc  is  not  to  be  broken;  and  all  other  /i^  and  are  0. 
Thus,  we  see  that  a  cut  is  identified  and  arcs  are  broken  in  that  cut  so  as  to  leave  as 
little  remaining  capacity  as  possible. 
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3.  Model  2  -  Network  with  Undirected  Arcs 


If  the  network  G’  =  (  N,  A’)  is  undirected,  we  would  have  to  convert  G’  into 
an  appropriate  directed  graph  in  order  to  use  Model  1.  In  the  maximum  flow 
problem,  no  flow  need  enter  the  source  nodes  and  no  flow  need  leave  the  sink  nodes 
so  that  any  arc  incident  to  a  source  node  can  be  represented  by  a  single  arc  directed 
out  from  the  source  node  and  any  arc  incident  to  a  sink  node  can  be  represented  by 
a  single  arc  directed  into  the  sink  node.  Other  undirected  arcs  in  G’  can  be  replaced 
by  two  independent  directed  arcs  in  anti-parallel,  each  directed  arc  with  the  capacity 
Ujj  and  resource  e,j  taken  from  the  original  network  G.  If  one  of  the  two  directed 
arcs  in  the  anti-parallel  pair  is  destroyed,  there  is  no  possibility  of  the  other  arc  being 
interdicted  simultaneously  requiring  a  total  of  le^j  units  of  resource  to  destroy  the 
undirected  arc.  Specifically,  if  one  directed  arc  (/,;)  is  destroyed,  we  have  a,  =  0  and 
Qj  =  1.  This  means  that,  for  the  other  directed  arc  {j,i),  the  term  {h^j  +  ^,7,^7}  will 
be  zero  without  violating  the  second  constraint  of  (2.5). 

Model  1  can  solve  the  problem  of  undirected  graph  using  the  above 
transformation,  but  introduces  a  large  number  of  variables  as  it  does  not  take 
advantage  of  the  fact  that  each  pair  of  directed  arcs  share  the  same  values  of 
capacity  and  destruction  effort.  Denote  an  undirected  arc  as  an  unordered  pair  of 
nodes  [/,;] 
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where  ij  e  G\  The  modified  model  becomes: 


Minimize  E  u^  h^- 
('VI  ■' 


(2.6) 


Subject  to 

0:  -  Or 
t  s 

=  1 

a.  -  a.  ^ 

N  *  ^T,;, 

>  0 

V  [/,/] 

a.  -  a.  + 

h,j  * 

>  0 

V  [/,;] 

<E, 

V  I 

>  0 

V  [iJ] 

f  (0,1) 

V  [iJ]  ,  / 

0. 

unrestricted 

V  / 
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C.  MINIMIZE  INTERDICTION  EFFORT 


1.  Problem  Statement 

The  attacker  has  sufficient  resources  E  to  stop  completely  the  flow  of 
material  in  the  supply  network.  The  time  taken  for  force  projection  to  arc  (/,/)  and 
recovery  of  resource  type  I  back  to  base  depends  on  a  factor  The  aim  of  the 
interdictor  is  to  destroy  the  set  of  arcs  so  as  to  cut  the  network  into  two  disjoint 
subsets  utilizing  minimum  resource  and  assigning  the  missions  in  such  a  way  which 
minimizes  the  maximum  mission  turnaround  time. 

2.  Model  3  -  Network  with  Directed  Arcs 

When  the  interdictor  has  only  one  type  of  resource,  the  minimum  resource 
required  to  reducing  the  network  flow  to  zero  and  the  associated  cut-set  can  be 
determined  by  solving  a  maximum  flow  problem  with  the  arc  capacities  equal  to  the 
interdiction  efforts.  The  minimal  capacity  cut-set  of  the  maximum  flow  problem 
yields  the  results  of  the  total  minimum  resource  and  the  set  of  arcs  to  be  interdicted. 

To  solve  the  problem  of  multiple  resources,  we  need  to  introduce  the  following 
modifications  to  Model  1: 

a.  Delete  all  terms  with  as  there  is  now  zero  network  flow, 

b.  Introduce  a  conversion  factor  w^j  for  each  resource  type  /;  the  conversion 

factor  being  the  measure  of  the  mission  turnaround  time  and  of  how  the  strike 
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planners  value  the  opportunity  cost  of  a  particular  resource  availability.  The 
formulation  becomes: 


Minimize  E  E  h> 

1  ('•;■) 

ijl 

Subject  to  Of,  - 

=  1 

>  0 

V  (/./) 

('•;) 

<E, 

V  / 

f  (0,1) 

V  (/,/)  ,  / 

unrestricted 

V  / 

(2.7) 
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3.  Model  4  -  Network  with  Undirected  Arcs 

For  an  undirected  network,  let  an  undirected  arc  be  represented  by  a  pair  of 
nodes  [/,;],  the  formulation  is  as  follows: 


Minimize  E  E 

1  l'.;l 

Subject  to  0,  -  ttj 

=  1 

a,  -  Q-^  +  E  y 

1 

>  0 

V  [ij] 

/ 

>  0 

V  [iJ] 

('w/1 

<E, 

V  / 

€  (0,1) 

V  [iJ] 

0. 

unrestricted 

V  / 

(2.8) 
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D.  MAXIMIZE  SHORTFALL  IN  PRIORITIZED  DEMAND 


1.  Problem  Statement 

The  enemy  wants  to  push  through  a  road  network  sufficient  flow  of 
materials  from  a  set  of  depots  to  meet  the  demands  of  the  field  units,  especially 
of  those  tactically  crucial  units  which  are  in  critical  shortage  of  supplies.  The 
problem  of  the  attacker  is  to  select  a  set  of  arcs  with  a  fixed  resource  E  in  order  to 
maximize  the  prioritized  shortfall  of  the  demands  measured  according  to  tactical 
values  r,  of  the  units. 

2.  Model  5  -  Network  with  Directed  Arcs 

In  this  model  the  additional  notation  used  is: 

.  the  set  of  supply  nodes, 

Ij  :  the  set  of  demand  nodes, 

4  :  the  set  of  transhipment  nodes. 

The  additional  data  required  is: 

5,  :  the  amount  of  supply  at  node  /. 

d,  :  the  amount  of  demand  at  node  /, 

r,  :  the  tactical  values  of  units;  the  higher  the  value  of  the  unit  the 

greater  is  its  tactical  importance. 

The  additional  variables  introduced  here  are: 

y*  :  the  amount  of  unused  supply  at  supply  node  /, 
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yf  :  the  amount  of  shortfall  at  demand  node  i. 
The  Max-Min  problem  formulation  is  as  follows: 


Maximize  Minimize  E  r-  y[ 
Y  £  r  it  I. 


(2.9) 


Subject  to 


= 

V  (/•  c  /,) 

E  Xy  -  E  x^, 

j  J 

=  0 

V  {i  f  /.) 

^^ij  -  -yi 

i  j 

=  -d- 
1 

V  {i  f  4) 

^ij 

^  »^,>(1  -  7yv) 

V  (/,/) ,  / 

>  0 

V  (»j) 

* 

y, 

>  0 

V  (/  £  /,) 

yi 

>  0 

V  (/  £  I,) 

where  X  =  {7,y,  1  7^,  «  (0,1)  V  (/J)  ,  I  ,and  E  7^/  <  £,  V  /  } 
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Let  a,  be  the  dual  variable  associated  with  the  flow  conservation  constraint  and 
the  capacity  constraint.  Converting  the  minimization  LP  into  its  dual,  we  have: 


Max  Max  E  5,0.  -  E  +  E  E  u..(l  -  7./)  (2.10) 

Y  £  r  I  £  /,  I  £  /^  /  (i.j) 


Subject  to  a.  -  a  +  E^^./  <0  V  (/,/) 

a,  <  0  V  (/  £  /^) 

-a.  <  r.  V  (/•  €  4) 

4^/  <  0  V  (/,/)  ,  / 


a.  unrestricted  V  /  €  (/^  ,  4) 

r  =  {  4/  I  4/  ^  (0^1)  V  (ij)  >  I  ,  ^  4,7y,  <  £/  V  /  } 


The  model  can  be  linearized  by  replacing  '^ith  h^j,  and  adding  the 

constraint  /jy,  <  4/  *  where  is  a  lower  bound  on  P^j,.  Looking  back  at  (2.8), 
since  P,^,  is  the  dual  variable  for  constraint  x^-  <  u-j{  1  -7^,)  ,  it  can  be  interpreted 

as  the  rate  of  change  of  the  objective  function  value  if  an  arc  is  interdicted.  Since 
is  negative  and  the  maximum  rate  of  change  in  the  objective  function  value  of  (2.8) 
is  max  r,  ,  it  follows  that  B  >  -max  r,  . 
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The  formulation  for  the  MIP  is  thus; 


Maximize 

E  5, or,  -  E  +  E  E 

It/,  1 « /^  /  (ij) 

(2.11) 

Subject  to 

a,  -  <  0 

V  (/,/) 

V  (/,/) ,  / 

(‘J) 

V  / 

y,,  f  (0.1) 

V  (/,;•)  ,  / 

o 

VI 

M  (/  £  /J 

-tt,  <  r, 

V  (/  £  /,) 

a.  unrestricted 

V  /€(/,./,) 

o 

VI 

V  (/,/),/ 

h,,  <  0 

V  (/J),/ 
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3.  Model  6  -  Network  with  Undirected  Arcs 

For  undirected  graph,  we  can  apply  a  transformation  similar  to  the  previous 
sections  to  the  network,  and  the  formulation  is  as  follows: 


Maximize 


Subject  to 


i  «  /, 


I  I'  Jl 


(2.12) 


a,  -  -  E  <  0 

<  0 


V  UJ] 
^  I  Ji 


I'd 


,yip  ^  0 

V  [/,/] ,  / 

VI 

\/  / 

7^v  f  (0.1) 

V  [/,/] ,  1 

a,  <  0 

V  [/  c  /J 

-a,  <  r, 

V  [/  f  /.] 

Oj  unrestricted 

V  /€[/.,/.] 

lA 

O 

V  [/,;•],/ 

o 

VI 

V  [/,;■]  ,  / 
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III.  DETERMINISTIC  MODELS  -  MULTICOMMODITY  NETWORKS 


A  military  logistics  supply  problem  usually  involves  the  expeditious  shipment 
of  certain  distinct  commodities,  such  as  ammunition,  fuel,  spare  parts,  food  and  etc, 
from  the  supply  depots  and  ammo  dumps  to  their  respective  destinations  along  the 
arcs  of  an  underlying  transportation  network.  This  scenario  results  in  the  well- 
known  multicommodity  flow  problem.  It  also  occurs  in  communication  systems, 
urban  traffic  systems  and  railway  problems,  as  well  as  in  many  others.  The 
commodities  interact  when  flowing  on  the  same  arc  either  by  competing  for  arc 
capacity,  or  by  causing  congestion. 

A.  MINIMIZE  MAXIMUM  FLOW 

1.  Problem  Statement 

In  this  operational  setting,  the  problem  of  the  strategic  strike  planners  can  be 
stated  as  follows: 

The  enemy  has  a  policy  of  forcing  through  the  unbroken  arcs  of  the  network 
the  flow  of  a  range  of  commodities.  Commodity  k  flows  from  source  node  5^ 
to  a  sink  node  The  aim  of  the  interdictor  is  to  reduce  the  sum  of  the 
maximum  flow  of  the  multiple  commodities  utilizing  a  fixed  amount  of 
resources  and,  if  given  sufficient  resources,  to  stem  the  flow  of  materials  with 
minimum  resource  and  mission  turnaround  time. 
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2.  Model  7  -  Network  with  Directed  Arcs. 

For  the  multicommodity  network  model,  the  following  notation  is  used: 
i,j  :  the  index  sets  for  the  nodes  of  the  network, 

k  :  the  index  set  lor  the  commodity, 

/  ;  the  index  set  for  the  types  of  strike  resources. 

The  given  data  for  the  model  are: 

u,j  :  the  combined  capacity  per  unit  time  in  a  standard  unit  for  the  arc 

(U)- 

The  variables  in  the  models  are: 

:  the  amount  of  flow  for  commodity  k  through  the  network, 

x,jk  :  the  amount  of  flow  for  commodity  k  in  arc  (/,/); 

:  the  binary  variable  which  is  1  if  arc  (/,/)  is  destroyed  by  strike 

resource  type  /,  and  is  otherwise  0. 
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rhe  formulation  of  the  interdiction  problem  from  the  attacker’s  viewpoint  is  as 


folio' vs: 


Minimize  Maximize  E 


Y  «  r 


Subject  to 


=1 
i 


h 

0 

-A 


V 


^.,k  ^  0 


if  i  =  s, 
if  /  f 
if  /  =  t, 

V  (/,/)  ,  / 

V  (/,/)  ,  k 


V  / 


r  =  {  7^/  I  7y/  €  (0,1)  V  (/,/)  ,  /  ,  and  E  e..,7yv  <  £/,  V  /  } 

(ij) 


(3.1) 


Let  be  the  dual  variable  associated  with  the  flow  conservation  equations,  and 
let  be  the  dual  variable  associated  with  the  combined  capacity  constraints.  By 


taking  the  dual  of  the  maximization  problem  the  Min-Max  problem  becomes: 


Minimize 
Y  c  r 

Minimize  E 
/ 

E  u  -  -iij,) 

('J) 

Subject  to 

=  1 

V  k 

+  E  $.j, 

>  0 

V  (/,;)  ,  k 

^iji 

>  0 

V  (/J)  ,  1 

unrestricted 

M  i  ,  k 

r  =  {7, 

,  1  hj,  ^  OJ) .  1 

,  and  E  e,j,i,j, 
('•» 

<  £,  V  /} 

The  above  model  can  be  linearized  by  replacing  -  7,^/)  with  /i,^,  and 
introducing  an  additional  constraint: 

^  (L/)  »  I 


where  is  a  known  upper  bound  on  and  a  suitable  value  is  =  1.  In  the 
multicommodity  flow  problem,  not  all  dual  variables  of  the  maximum  flow  problem 
need  to  be  0  or  1.  There  is  no  guarantee  that  the  optimum  solution  will  be  integer 
(unless  we  have  a  completely  planar  network,  Sakarovitch  1966). 
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The  full  model  is: 


Minimize  E  E  w,  /i,y 
/  Uj) 


Subject  to 


=  1 

V  k 

0,1, 

>  0 

V  (/,;)  ,  k 

^ijl 

>  0 

V  (/,/) ,  / 

Oj) 

<E, 

V  / 

c  (0,1) 

V  (/ J) ,  / 

hi 

>  0 

V  (/./) ,  1 

hi 

>  0 

V  (/J) .  / 

<^ik 

unrestricted 

M  i  ,k 

(3.3) 
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3.  Model  8  -  Network  with  Undirected  Arcs. 


In  the  undirected-arc  multicommodity  flow  problem,  the  material  can  flow  in 
either  direction  of  an  arc  so  that  the  total  flow  is  limited  by  the  sum  of  flows  of 
different  commodities  in  both  directions.  An  undirected  arc  in  this  type  of  network 
can  be  modeled  by  two  anti-parallel  directed  arcs  with  a  joint  capacity  constraint. 
Flows  of  same  commodity  in  opposite  direction  will  cancel,  while  flows  of  different 
commodities  are  cumulative  regardless  of  direction.  The  arc  capacity  constraint  for 
undirected  arc  [/,//  is  then: 


The  Min-Max  problem  from  the  attacker’s  viewpoint  is; 


Minimize 
Y  £  r 


r  ={y.j, 


Maximize  E /j,.  (3.4) 

k 

Subject  to 


j  J 


fk 

0 

-fk 


if  i  =  s. 


if  i  s„t, 


if  /  = 


^  (  -V-  ^  ^j.k  )  ^ 


Xijk  ^  0 


V  k 

V  [i,j]  ,  I 

V  [/,;]  ,  k 


I  7  €  (0,1)  V  [ij]  ,  /  ,  and  E  e../7yv  <  E,,M  I  } 

I'jl 
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Let  an  undirected  arc  be  denoted  by  [i,j].  Following  the  same  methodology  for 
the  development  of  Model  6,  the  formulation  for  the  undirected-arc  multicommodity 


network  problem  becomes: 


Minimize 

E  E 

1  liJ] 

Subject  to 

-^sk  *  ^tk 

=  1 

V  k 

- 

>  0 

V  [/,/]  ,  k 

-  O/t  * 

1 

>  0 

V  [/,/]  ,  k 

>  0 

V  [/,;•]  ,  1 

[‘■J] 

<E, 

V  / 

€  (0,1) 

V  [/J]  ,  1 

>  0 

V  [/,;] ,  / 

fhj, 

>  0 

V  [/,/],/ 

^ik 

unrestricted 

V  /  ,  k 

(3.5) 
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B.  MINIMIZE  INTERDICTION  EFFORT 


1.  Problem  Statement 

The  attacker  has  sufficient  resources  to  completely  stop  the  flow  of  multiple 
commodities  through  the  supply  network.  The  aim  of  the  attacker  is  to  create  with 
the  minimum  resource  and  minimum  mission  turnaround  time  a  discoimecting  set 
which  is  defined  as  a  set  of  arcs  whose  removal  from  the  network  destroys  all  paths 
from  each  source  to  its  respective  sink  . 

2.  Model  9  -  Network  with  Directed  Arcs 

To  solve  the  above  problem,  we  need  to  introduce  the  following  modifications 
to  Model  6  in  (3.3): 

a.  since  there  are  now  sufficient  resources  and  we  assume  that  we  have 
solved  Model  6,  all  /i,y/  terms  in  the  objective  function  will  be  zero  indicating 
zero  flow  in  the  network.  All  /i.y,  terms  may  then  be  deleted  from  the  model, 

b.  the  objective  function  of  the  new  model  should  now'  be  the  amount  of 
total  resources  weighted  by  the  factor  mv,  which  is  the  measure  of  the 
turnaround  time  for  mission  of  type  /  on  arc  (/,;). 
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The  model  is  as  follows: 


Minimize  E  E  'i  , 

1  oj) 

Subject  to 

V  k 

,  -  Q ,  +  E^  ,  >  0 

ik  ik  ^ijl  - 

1 

V  (/,;■)  .  ^ 

-p.ji  *  ^  0 

V  (/./)  ,  / 

Oj) 

\/  / 

?  0 

V  (/J)  ,  / 

1,,  ‘  (0.1) 

V  (/,/) ,  / 

unrestricted 

M  i  ,  k 

(3.6) 
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3.  Model  10  -  Network  with  Undirected  Arcs 

For  the  undirected  graph,  the  formulation  becomes: 


Minimize  E 
1 

Subject  to 

V  ‘'jl 
['Jl 

=  1 

V  k 

>  0 

V  [/.;•]  ,  k 

>  0 

V  (/,;■]  ,k 

>  0 

V  [/,/]  ,  / 

<-E, 

V  / 

^,ji 

>  0 

V  It  Jl  ,  1 

f  (0,1) 

V  [/,/] .  1 

unrestricted 

i  ,  k 

IV.  PROBABILISTIC  MODELS 

In  this  chapter,  we  will  present  probabilistic  models  for  network  interdiction. 
The  models  to  be  presented  are  applicable  to  infiltration  and  counterinsurgency 
situation. 

A.  ONE  EVADER  VERSUS  ONE  INTERDICTOR 

1.  Problem  Statement 

An  evader  attempts  to  travel  from  the  source  node  5  of  a  supply  network 
to  the  sink  node  t.  An  interdictor  would  like  to  stop  his  opponent  by  positioning  with 
probability  ii,j  on  arc  (/,/)  so  that  the  probability  of  successful  infiltration  by  the 
evader  is  minimized.  The  probability  that  the  interdictor  detects  the  evader  given 
that  the  interdictor  positions  himself  on  arc  (/,/)  and  the  evader  traverses  on  arc  (/,;) 
is  pij.  (If  detection  is  certain,  p^j  could,  instead,  represent  a  reward  for  catching  the 
evader  on  arc  (/,/).)  Knowing  the  interdictor’s  strategy,  the  evader  will  have  to 
develop  an  optimum  strategy  by  selecting  path  k  with  probability  6^  so  as  to  minimize 
the  probability  of  interception. 

2.  Model 

For  this  model,  the  indices  used  are: 

i,j  :  nodes. 

The  data  of  the  model  are: 
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Pij  :  probability  of  successful  interception  given  that  the  evader  has  chosen 
to  cross  the  arc  and  the  interdictor  has  placed  himself  there. 

The  decision  variables  are: 

K,j  :  the  probability  of  the  interdictor  deploying  at  arc  (/,;), 

(f),j  :  the  probability  of  the  evader  traversing  arc  (;j). 

The  Min-Max  formulation  is  as  follows: 


Minimize  Maximize  E 

♦  It  (i.j) 


(4.1) 


Subject  To 


E  X  =  1 

UJ) 


E  0..  -  E  0^,.  = 

j  j 


1 

0 

-1 


>  0 


if  i  =  s 
if  i  f  s  OT  t 
if  i  =t 

V  (ij) 

V  (/,;) 


The  aim  of  the  interdictor  is  to  maximize  the  total  reward  while  the  aim  of  the 
evader  is  to  minimize  the  total  reward.  The  first  constraint  of  (4.1)  ensures  that  the 
probability  of  deployment  by  the  interdictor  is  one.  The  second  set  of  constraints 
ensure  that  only  one  evader  is  traversing  the  network.  Let  v  be  the  dual  variable 
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associated  with  the  first  constraint  of  (4.1)  and  taking  the  dual  of  the  maximization 
problem,  we  have: 


Minimize 


(4.2) 


Subject  To 


V  -  >  0 


E  -  E  (t)j.  = 
J  j 


K.J  >  0 


V  (i,j) 

1  if  /  =  5 

0  if  ^  5  or  r 

-1  if  /  =  t 

V  (/,/) 

V  (ij) 

V  unrestricted 


<!>,>- Q 


The  optimal  solution  will  occur  for  that  value  of  v  such  that  the  minimum  cut 
capacity  (  with  arc  capacities  at  v/pjj  )  is  1.  The  dual  prices  of  the  first  set  of 
constraints  in  the  optimum  solution  will  be  the  optimum  values  of  and  are  positive 
on  the  minimum  cut-set.  The  optimal  strategy  of  the  evader  can  be  recovered  as 
follows: 

(1)  Let 

if  path  k  contains  arc  (/,;) 
otherwise 


a 


0 
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and  6*  be  the  probability  that  the  evader  chooses  path  k. 


(2)  Having  known  the  optimal  values  of  by  solving  (4.2),  can  be 
determined  by  solving  the  following  set  of  equations: 

^  =  <l^ij  ^  (iJ) 

k 


As  Lawler  1976  points  out,  there  exists  a  solution  of  such  a  system  with  at  most 
l£l  positive  6*.  In  fact,  a  set  of  l£l  or  fewer  arcs  which  will  satisfy  this  system  can 
be  derived  as  follows:  Start  with  the  optimal  "flow  of  probability"  in  the  system  from 
(4.2).  Discard  any  arcs  with  zero  flow.  Find  a  directed  path  kj  from  s  to  t  and 
suppose  the  minimum  flow  arc  on  that  arc  has  flow  Xj.  Then,  associated  with  this 
path  k;,  6^^  =  x^  .  Subtract  X;  units  of  flow  from  arcs  along  the  path  and  delete  any 

arcs  with  no  flow  remaining.  Repeat  this  process  at  most  l£l  -  1  more  times  until 
all  flow  has  been  allocated  to  at  most  l£l  paths. 
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B.  ONE  EVADER  NTRSUS  MULTIPLE  INTERDICTORS 


1.  Problem  Statement 

A  team  of  n  interdictors  is  deployed  in  the  network  in  an  attempt  to  stop  an 
evader  from  traversing  from  the  source  node  to  the  sink  node.  Tiie  evader  has  the 
probability  0^  of  choosing  to  travel  on  arc  At  arc  (ij),  at  most  one  interceptor 
/  may  be  deployed,  with  the  probability  of  tty,.  The  reward  to  an  interdictor  of 
deploying  at  arc  (/,/)  given  that  the  evader  has  chosen  to  traverse  that  arc  is  The 
objective  of  the  interdictor  team  is  to  maximize  the  total  reward  while  that  of  the 
evader  is  to  minimize  it. 


2.  Model 

In  this  model,  the  additional  index  used  is: 

/  :  interdictor  number  which  is  l...n. 

The  additional  decision  variables  are: 

Tty,  :  the  probability  of  interdictor  /  deploying  at  arc  (/,/), 

S,ji  :  the  binary  variable  which  is  1  if  interdictor  I  is  deployed  at  arc  (/,;) 


and  is  otherwise  0. 


The  Max-Min  formulation  is  as  follows; 


Maximize  Minimize  E  p ^6 

«  .  S  <*>  iij)  ’  i  ’ 


Subject  To 


E  -  E  (t>j,  = 
J  J 


1 

0 

-1 


E  <  1 
/ 


E  IT  <  1 

UJ) 


^ij,  ^  0 


S,j,  e  (0,1) 


(4.3) 


if  1=5 
if  i  f  s  OT  t 
if  i  =  t 

V  (/,/)  ,  / 

V  (/,/) 

V  / 

V  (/,/) ,  / 

V  (/,/) 

V  UJ) ,  / 


The  objective  of  the  interdictor  is  to  manipulate  Tr^,  and  so  that  the  total 
reward  would  be  maximized.  The  evader  aims  to  minimize  this  reward,  manipulating 
(j),j.  The  objective  function  is  correct  since,  at  any  one  time,  only  one  interdictor 
would  be  present  at  an  arc.  This  is  so  due  to  the  second  and  third  constraints  of 
(4.3).  This  formulat’on  does  restrict  the  interdictors  in  that  if  interdictor  /  ever 
interdicts  arc  (/,/),  no  other  interdictor  1=1  can  ever  interdict  that  arc. 
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Taking  the  dual  of  the  minimization  LP  of  the  above  Max-Min  problem  and 
letting  a,  be  the  dual  variable  associated  with  the  first  constraint  of  (4.3),  we  have: 


Maximize 
Subject  To 


-  a,  (4.4) 


<  0 

V  i 

• 

- «,)/ 

<  0 

V  (/,/) ,  / 

♦ 

1  ■' 

<  1 

V  (/,;•) 

E  TT 

(ij) 

<  1 

V  / 

>  0 

V  (/,/) ,  / 

>  0 

V  ('J) 

£  (0,1) 

V  (/,;)  ,  / 
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V.  COMPUTATIONAL  RESULTS 

In  this  chapter,  a  network  interdictor’s  problem  which  originally  motivated  the 
writing  of  this  thesis  is  solved.  GAMS  (General  Algebraic  Modeling  System) 
programs  are  developed  and  lun  on  an  Amdahl  5990  mainframe  computer.  The 
integer  programs  are  solved  using  ZOOM  (Zero/One  Optimization  Methods)  which 
is  available  in  GAMS.  The  network  to  be  solved  is  an  undirected  network,  and  the 
relevant  models  tested  are: 

a.  Model  2  -  minimize  the  maximum  flow  (single  commodity  ), 

b.  Model  4  -  minimize  the  interdiction  effort  (single  commodity), 

c.  Model  8  -  minimize  the  maximum  flow  (multicommodity), 

d.  Model  10  -  minimize  the  interdiction  effort  (multicommodity). 

A.  SCENARIO 

The  road  transportation  network  in  question  is  as  shown  in  Figure  1.  It  is  a 
land  transportation  network  system  found  in  Southeast  Asia  and  consists  of  1 12  nodes 
and  176  undirected  arcs.  The  capacities  of  the  arcs  in  term  standard  truck -loads 
per  day  as  shown  in  Table  1  which  is  in  Appendix  A.  The  land  transportation 
network  in  question  presents  itself  as  a  realistic  scenario  for  a  strategic  strike  planner 
to  solve.  The  locations  of  enemy  front-line  units,  its  supply  depots  and  the 
orientation  of  the  network  are  as  shown  in  Figure  1. 
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1.  Opposing  Forces 

The  network  interdictor  has  at  his  disposal  the  following  types  of  assets: 


a.  Ground  Attack  Fighter  Squadrons.  Fighter  aircraft  laden  with 
weapons  such  as  laser-guided  bombs  and  other  precision  guided  munitions 
which  are  suitable  for  destructions  of  bridges. 

b.  Naval  Gunfire  Support.  Naval  combatants  can  be  deployed  along  the 
coasts  to  provide  gun  fire  support. 

c.  Special  Forces.  This  consists  of  highly  specialized  and  trained 
personnel  operating  autonomously  in  small  teams  deep  in  enemy’s 
territories  to  carry  out  clandestine  sabotage  missions. 

2.  Weapons  Efforts 

Against  a  particular  target,  the  suitability  of  a  weapon  system  and  the  required 
efforts  depend  on: 

a.  the  types  of  defenses  which  the  attacker  is  expected  to  encounter, 

b.  the  effective  range  of  the  weapon, 

c.  in  the  case  of  aerial  attacks,  the  degree  of  difficulty  in 
detecting  and  identifying  targets,  and, 

d.  the  size  of  target. 

The  resources  required  to  destroy  each  target  is  given  in  Table  1.  The  data 
presented  represent  a  realistic  range  of  possibilities  but  are  hypothetical  and  only 
meant  for  the  purpose  of  testing  the  model.  Accurate  resource  data  could  be 
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obtained  from  detailed  mission  planning  for  each  target  with  the  aid  of  the  Joint 
Munitions  Effectiveness  Manual  (JMEM),  1991. 

B.  MINIMUM  CUT-SET  -  SINGLE  COMMODITY  NTITVVORK 

In  this  scenario,  the  enemy  wishes  to  supply  his  front-line  troops  at  nodes  [108] 
through  [112]  by  forcing  the  flow  of  a  single  commodity  from  the  supply  depots  at 
nodes  [1]  through  [7]  through  the  network.  The  maximum  flow  of  the  un-interdicted 
network  is  190  units  and  the  minimum  cut-set  is  as  shown  in  Figure  2. 

C.  MODEL  2  -  MINIMIZING  MAXIMUM  FLOW  -  SINGLE  COMMODITY 
NETWORK 

For  this  model,  the  GAMS  program  that  solves  the  problem  is  in  the  Appendix. 
Given  that  the  interdictor  has  4  aircraft,  4  naval  gunfire  support  (NGFS)  units  and 
4  teams  of  special  forces  with  which  to  interdict  the  network,  the  maximum  flow  is 
reduced  to  20  units.  The  arcs  to  be  interdicted  and  the  force  allocation  are  as 
follows; 


Arcs 

Weapons  Allocation 

1. 

(84,85) 

2  Aircraft 

2. 

(83,86) 

4  NGFS 

3. 

(54,87) 

2  Aircraft 
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Arcs 

Weapons  Allocation 

4. 

(53,52) 

2  Special  Forces 

5. 

(51,94) 

2  Special  forces 

The  cut-set  is  as  shown  in  Figure  3.  In  this  cut-set,  arc  (50,52)  is  the  only  arc 
not  interdicted,  and  allows  a  flow  of  20  units.  The  model  contains  356  single 
equations,  817  real  variables  and  528  discrete  variables.  The  generation  time  and 
execution  times  are  1.58  seconds  and  1.66  seconds  respectively.  The  generation  time 
is  the  time  spent  preparing  the  model  for  solution,  while  the  execution  time  is  the 
time  used  after  the  syntactic  check  is  finished,  including  the  time  spent  generating  the 
model. 
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D.  MODEL  4  -  MINIMIZE  INTERDICTION  EFFORT  -  SINGLE  COMMODITY 


NETWORK 

Given  that  the  interdictor  has  sufficient  resources  to  stop  completely  the  flow 
of  material  in  the  supply  network,  and,  for  simplicity,  assuming  equal  turnaround 
times,  the  minimum  effort  required  and  the  force  allocation  are  as  follows: 


Arcs 

Weapons  Allocation 

1. 

(84.85) 

2  Aircraft 

2. 

(83,86) 

2  Aircraft 

3. 

(54,87) 

2  Aircraft 

4. 

(52,82) 

4  Special  Forces 

5. 

(51,94) 

2  Special  Forces 

Figure  4  shows  the  location  of  strikes  and  the  force  allocation.  The  model  has  356 
equations,  641  real  variables  and  628  discrete  variables.  It  takes  1.59  seconds  to 
generate  and  1.66  seconds  to  execute  the  model. 
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E.  MODEL  8  -  MINIMIZE  NUXIMLM  FLOW  -  MULTICOMMODITY 


NETWORK 

For  the  multicommodity  flow  scenario,  we  assume  that  the  enemy  wishes  to 
supply: 

1.  the  unit  at  node  [112]  with  w-ater  from  the  water  "source"  at  node  [90], 

2.  the  unit  at  node  [108]  with  ammunition  from  the  ammunition  dump  at 

node  [6], 

3.  the  unit  at  node  [112]  with  fuel  from  the  fuel  depot  at  node  [5]. 

The  maximum  multicommodity  flow  is  260  units  which  is  the  sum  of  flow  of  the 
three  commodities  in  terms  of  standard  truckloads  per  day.  With  the  available  assets 
of  4  aircraft,  4  NGFS  units  and  4  teams  of  special  forces,  the  maximum  flow  is 
reduced  to  100  units.  Figure  5  shows  the  locations  of  attack  and  the  following  force 
allocation: 


Arcs 

Weapons  Allocation 

1. 

(51,94) 

2  Special  Forces 

2. 

(54,55) 

2  Aircraft 

3. 

(54,83) 

2  Aircraft 

4. 

(83,86) 

4  NGFS 

5. 

(84,85) 

2  Special  Forces 
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Note  that  the  unit  at  node  [112]  will  receive  80  units  of  water  supply  and  20  units  of 
fuel  and  ammunition  will  flow  through  the  network.  The  model  has  1588  equations, 
1921  real  variables  and  528  discrete  variables.  It  takes  6.01  seconds  to  generate  and 
6.11  second  to  execute  the  model. 


F.  MODEL  10  -  MINIMIZE  INTERDICTION  EFFORT  -  MULTICOMMODITY 
NETWORK 

To  disconnect  all  the  units  concerned  from  their  respective  sources  of  supply, 
the  minimum  interdiction  effort  is  14  units  and  the  force  allocation  is  as  follows: 


Arcs 

Weapons  Allocation 

1. 

(5,12) 

4  Aircraft 

2. 

(6,9) 

2  Aircraft 

3. 

(97,112) 

8  NGFS 

As  shown  in  Figure  6,  the  interdiction  efforts  are  channelled  to  the  arcs  incident  to 
the  fuel  depot,  ammunition  dumps  and  the  tactical  unit  requiring  water  supply.  If  the 
enemy  augments  the  defenses  of  these  arcs  by  increasing  their  destruction  efforts  e,^, 
to  an  extent  that  it  is  no  longer  possible  for  the  attacker  to  strike  these  arcs  ( in  this 
example,  we  make  >  £,  for  these  three  arcs),  we  will  see  a  distinct  shift  in  the 
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attacker’s  strategy.  Figure  7  shows  the  new  strategy  adopted  by  the  attackers.  The 
reallocation  of  the  strike  resources  is  as  follows: 


Arcs 

Weapons  Allocation 

1. 

(84,85) 

2  Aircraft 

2. 

(83,86) 

2  Aircraft 

3. 

(54.87) 

2  Aircraft 

4. 

(52,94) 

4  Special  Forces 

5. 

(51,94) 

2  Special  Forces 

6. 

(90,89) 

4  NGFS 

7. 

(90,91) 

4  Aircraft 

The 

variables. 


model  has  1588  single  equations,  1393  real  variables  and  528  discrete 
It  takes  5.72  seconds  to  generate  and  5.82  seconds  to  execute  the  model. 
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VI.  SUMMARY  AND  EXTENSIONS 


A.  SUMMARY 

In  this  report,  we  present  two  main  types  of  model,  vis.,  the  deterministic  and 
probabilistic  models,  for  the  analysis  of  strategic  strike  against  a  land  communication 
network. 

In  the  deterministic  models,  the  technique  of  mixed  integer  programming  is 
used  and  has  been  shown  to  be  quite  versatile  in  that  generalizations  to  include, 
among  other  features,  multiple  strike  resources,  are  easy.  Under  this  category,  the 
main  models  we  have  developed  are: 

a.  Single  Commodity  Network.  The  aim  of  the  network  interdictor  is  to 
minimize  the  maximum  throughput  of  material  supply  by  allocating  a  fbced 
amount  of  multiple  resources  to  interdict  the  most  "profitable"  arcs  of  the 
underlying  transportation  network.  If  there  are  sufficient  resources  to  stop 
completely  the  flow  of  materials,  the  problem  becomes  one  of  allocating,  in  an 
optimal  manner,  the  minimum  effort  over  the  multiple  strike  resources. 
Another  way  of  degrading  the  performance  of  the  supply  network  is  to 
maximize  the  prioritized  shortfalls  of  the  demands  at  the  tactical  units. 

b.  Multicommodity  Network.  In  this  situation,  the  enemy  wishes  to  send  a 
few  types  of  commodities  from  their  supply  depots  to  the  respective  tactical 
units  requiring  specialized  support.  The  attacker’s  problem  is  either  to 
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minimize  the  sum  of  the  maximum  flow  of  the  multiple  commodities  or  to 
utilize  minimum  resources  to  achieve  a  disconnecting  set,  which  severs  the 
paths  connecting  all  sources  to  their  respective  sinks. 

Another  type  of  model  presented  is  a  probabilistic  model,  which  is  applicable 
to  the  interdiction  of  a  lightly  travelled  network  in  which  the  arc  capacity  is  not  a 
factor.  The  objective  of  a  single  interdictor  or  a  team  of  interdictors  is  to  minimize 
the  probability  of  infiltration  by  a  single  evader  through  a  network  while  the 
objective  of  the  evader  is  just  the  opposite. 

B.  EXTENSIONS 

There  are  a  few  useful  model  extensions  which  warrant  further  research. 

1.  Deterministic  Models 

a.  Partial  Arc  Destruction.  One  possible  model  extension  would  be  to  allow 
for  partial  arc  destruction,  by  assuming  a  linear  or  piece-wise  linear  damage 
function. 

b.  Co-ordinated  Strikes.  Certain  missions  may  require  the  co-ordinated 
execution  of  a  few  upes  of  strike  resources.  Examples  include  the  employment 
of  ground  laser  designator  to  illuminate  targets  for  fighter  bombers,  and  target 
spotting  for  naval  gun  fire  support.  The  model  presented  can  be  modified  to 
incorporate  this  feature. 
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2.  Probabilistic  Models 


a.  Multiple  Evaders  Versus  Multiple  Interceptors.  It  will  be  useful  to  extend 
the  model  to  include  the  situation  where  multiple  interceptors  are  available  to 
defend  the  network  against  multiple  evaders. 

b.  Resource  Constraint.  The  issue  of  resource  constraint  over  time  can  be 
addressed  to  improve  model  realism. 
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APPENDIX  A.  NETWORK  DATA 


Arc 

Capacity 

Weapons  Efforts 

Aircraft 

NGFS 

Sp  Forces 

1. 

100 

12 

10 

15 

2. 

(2,20) 

40 

14 

8 

15 

3. 

(3,19) 

100 

14 

10 

15 

4. 

(4,14) 

no 

6 

12 

15 

5. 

(5,12) 

80 

4 

9 

15 

6. 

(6,9) 

130 

2 

15 

15 

7. 

(7,8) 

30 

12 

6 

15 

8. 

(8,43) 

50 

10 

6 

15 

9. 

(8,41) 

160 

14 

18 

15 

10. 

(8,9) 

140 

12 

15 

15 

11. 

(9,41) 

100 

16 

15 

15 

12. 

100 

18 

15 

15 

13. 

80 

12 

9 

15 

14. 

90 

16 

10 

15 

15. 

(11,37) 

100 

14 

12 

15 

16. 

(11,12) 

50 

16 

28 

15 

17. 

(12,36) 

50 

4 

29 

2 

18. 

(12,13) 

40 

6 

24 

4 

19. 

(13,33) 

80 

4 

28 

4 

20. 

(13,31) 

100 

8 

22 

5 

21. 

(13,14) 

20 

4 

23 

15 

22. 

(14,15) 

40 

4 

25 

15 

23. 

(15,31) 

50 

4 

28 

15 

24. 

(15,30) 

50 

2 

28 

15 

25. 

(15,28) 

60 

4 

6 

15 

26. 

(15,16) 

70 

6 

9 

2 

27. 

(16,23) 

130 

4 

14 

5 

28. 

(16,17) 

120 

6 

1 

15 

4 
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Capacity 


Weapons  Efforts 


Aircraft  NGFS 


Sp  Forces 


68 


Capacity 


Weapons  Efforts 


Aircraft  NGFS 


Sp  Forces 


Capacity 


Weapons  Efforts 


Aircraft  NGFS 


93. 

94. 

■ 

95. 

(57,59) 

30 

96. 

(57,58) 

60 

97. 

(58,59) 

50 

98. 

(59,61) 

150 

99. 

(59,62) 

140 

100. 

(60,61) 

20 

101. 

(62,63) 

70 

102. 

(62,64) 

90 

103. 

(63,71) 

40 

104. 

(63,64) 

30 

105. 

(64,70) 

80 

106. 

(64,65) 

80 

107. 

(65,70) 

50 

108. 

(65,67) 

40 

109. 

(65,66) 

10 

110. 

(66,67) 

120 

111. 

(67,68) 

130 

112. 

(67,69) 

60 

113. 

(69,76) 

70 

114. 

(69,70) 

40 

115. 

(70,71) 

60 

116. 

(71,76) 

70 

117. 

(71,72) 

80 

118. 

(72.75) 

40 

119. 

(72,73) 

140 

120. 

(73,74) 

150 

121. 

(74,79) 

120 

122. 

(74,76) 

60 

123. 

(74.75) 

30 

124. 

(75,76) 

20 

Sp  Forces 


5 

2 

2 

2 

2 

15 

15 

15 

5 

5 

5 

2 

5 

5 

2 

2 

2 

8 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 
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Arc 

Capacity 

Weapons  Efforts  || 

Aircraft 

NGFS 

Sp  Forces 

IB 

(76,78) 

5 

15 

IB 

(76,77) 

4 

15 

127. 

(77,78) 

80 

WS^M 

8 

15 

128. 

50 

rn^rn 

6 

15 

129. 

(78,79) 

30 

12 

5 

15 

130. 

(79,81) 

140 

14 

18 

15 

131. 

90 

12 

5 

132. 

30 

14 

4 

2 

133. 

40 

10 

4 

3 

134. 

(80.84) 

50 

12 

8 

2 

135. 

(81,83) 

50 

2 

8 

2 

136. 

(81,82) 

60 

2 

8 

5 

137. 

(83,86) 

30 

2 

4 

2 

138. 

(83,84) 

90 

2 

10 

6 

139. 

(84,85) 

80 

2 

8 

2 

140. 

(85,106) 

40 

4 

8 

15 

141. 

(85,86) 

70 

4 

8 

15 

142. 

(86,89) 

10 

8 

4 

15 

143. 

(86,87) 

160 

8 

20 

15 

144. 

(87,88) 

180 

6 

24 

15 

145. 

(88,89) 

40 

6 

4 

15 

146. 

(89,103) 

30 

4 

4 

15 

147. 

(89,100) 

90 

2 

13 

15 

148. 

(89,90) 

40 

6 

4 

15 

149. 

(90,91) 

60 

4 

5 

15 

150. 

(91,99) 

90 

18 

11 

15 

151. 

(91,98) 

60 

16 

5 

15 

152. 

(91,97) 

40 

14 

4 

15 

153. 

(91,96) 

140 

14 

12 

15 

154. 

(91,92) 

150 

10 

12 

15 

155. 

(92.93) 

150 

10 

12 

15 

156. 

(93.96) 

30 

14 

4 

15 
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Arc 

Capacity 

Weapons  Efforts 

Aircraft 

NGFS 

Sp  Forces  | 

12 

4 

158. 

12 

8 

159. 

12 

10 

15 

160. 

90 

12 

10 

15 

161. 

80 

14 

8 

15 

162. 

90 

12 

10 

15 

163. 

80 

12 

8 

15 

164. 

(99,100) 

60 

12 

6 

15 

165. 

(102,103) 

50 

12 

6 

15 

166. 

(103,107) 

90 

12 

10 

15 

167. 

(103,104) 

150 

10 

14 

15 

168. 

(104,107) 

120 

10 

12 

15 

169. 

(104,105) 

130 

16 

12 

16 

170. 

(104,106) 

50 

10 

8 

14 

171. 

(105,106) 

80 

10 

8 

14 

172. 

(106,108) 

100 

14 

10 

15 

173. 

(107,109) 

90 

14 

8 

14 

174. 

(101,110) 

100 

12 

10 

16 

175. 

(91,111) 

100 

12 

10 

16 

176. 

(97,112) 

80 

12 

8 

14 
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APPENDIX  B.  COMPUTER  PROGRAM 


GAMS  COMPUTER  LISTING 

08/04/91  10:31:48  PAGE  1 

2.19  IBM  CMS 


8 

9  OPTIONS 

10  LIMCOL  =  0,  UMROW  =  0,  SOLPRINT  =  OFF  ,  DECIMALS 

11  RESLIM  =10000,  ITERLIM  =  90000,  OPTCR  =  0.05; 

12 

13  * - DEFINITIONS  AND  DATA . . 

14 

15  SET 

16  I  nodes  in  the  network  /I*  112/ 

17  L  types  of  asset  /AIR,NGS,SF/; 

18 

19  ALIAS  (I,J): 

20 


39 

12  .36 

50 

40 

12  .13 

40 

41 

13  .33 

80 

42 

13  .31 

100 

43 

13  .14 

20 

44 

14  .15 

40 

45 

15  .31 

50 

46 

15  .30 

50 

47 

15  .28 

60 

48 

15  .16 

70 

49 

16  .23 

130 

50 

16  .17 

120 

51 

17  .23 

120 

52 

17  .22 

110 

53 

17  .18 

120 

54 

18  .19 

120 

55 

19  .22 

30 

56 

19  .20 

50 

57 

20  .21 

130 

58 

21  .22 

40 

59 

21  .24 

50 

60 

22  .23 

40 

08/04/91  10:31:48 
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61 

23  .28 

80 

62 

23  .26 

70 

63 

23  .24 

90 

64 

24  .25 

50 

65 

25  .26 

60 

66 

25  .27 

130 

67 

25  .68 

20 

68 

26  .27 

50 

69 

27  .28 

60 

70 

27  .62 

40 

71 

27  .66 

90 

72 

28  .29 

50 

73 

29  .61 

70 

74 

29  .30 

50 

75 

30  .61 

60 

74 


76 

30  .60 

70 

77 

30  .31 

100 

78 

31  .32 

120 

79 

32  .33 

150 

80 

33  .60 

140 

81 

33  .34 

60 

82 

34  .38 

20 

83 

34  .35 

140 

84 

35  .36 

160 

85 

35  .37 

30 

86 

38  .60 

60 

87 

38  .57 

40 

88 

38  .39 

20 

89 

39  .45 

80 

90 

39  .40 

90 

91 

41  .42 

40 

92 

42  .44 

90 

93 

42  .43 

60 

94 

43  .47 

30 

95 

44  .45 

140 

96 

44  .46 

160 

97 

46  .47 

180 

98 

47  .48 

130 

99 

48  .49 

140 

100 

49  .50 

30 

101 

49  .51 

150 

102 

51  .94 

50 

103 

50  .52 

20 

104 

52  .94 

90 

105 

52  .53 

30 

106 

53  .54 

50 

107 

54  .87 

30 

108 

54  .83 

70 

109 

54  .55 

30 

110 

55  .81 

150 

111 

55  .82 

140 

112 

55  .56 

30 

113 

56  .82 

70 

K4 

56  .73 

80 

115 

56  ,57 

90 

116 

57  .72 

40 

08/04/91  10:11:48  PAGE  3 
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GAMS 


117 

57  .59 

30 

118 

57  .58 

60 

119 

58  .59 

50 

120 

59  .61 

150 

121 

59  .62 

140 

122 

60  .61 

20 

123 

62  .63 

70 

124 

62  .64 

90 

125 

63  .71 

40 

126 

63  .64 

30 

127 

64  .70 

80 

128 

64  .65 

80 

129 

65  .70 

50 

130 

65  .67 

40 

131 

65  .66 

10 

132 

66  .67 

120 

133 

67  .68 

130 

134 

67  .69 

60 

135 

69  .76 

70 

136 

69  .70 

40 

137 

70  .71 

60 

138 

71  .76 

70 

139 

71  .72 

80 

140 

72  .75 

40 

14i 

72  .73 

140 

142 

73  .74 

150 

143 

74  .79 

120 

144 

74  .76 

60 

145 

74  .75 

30 

146 

75  .76 

20 

147 

76  .78 

50 

148 

76  .77 

30 

149 

77  .78 

80 

150 

77  .80 

50 

151 

78  .79 

30 

152 

79  .81 

140 

153 

79  .80 

90 

154 

80  .81 

30 

155 

8C  .83 

40 

76 


156 

80  .84 

50 

157 

81  .83 

50 

158 

81  .82 

60 

159 

83  .86 

30 

160 

83  .84 

90 

161 

84  .85 

80 

162 

85  .106 

40 

163 

85  .86 

70 

164 

86  .89 

10 

165 

86  .87 

160 

166 

87  .88 

180 

167 

88  .89 

40 

168 

89  .103 

30 

169 

89  .100 

90 

170 

89  .90 

40 

171 

90  .91 

60 

172 

91  .99 

90 
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173 

91  .98 

60 

174 

91  .97 

40 

175 

91  .96 

140 

176 

91  .92 

150 

177 

92  .93 

150 

178 

93  .96 

30 

179 

93  .94 

40 

180 

94  .95 

80 

181 

95  .97 

90 

182 

95  .96 

90 

183 

98  .99 

80 

184 

99  .101 

90 

185 

99  .102 

80 

186 

99  .100 

60 

187 

102.103 

50 

188 

103.107 

90 

189 

103.104 

150 

190 

104.107 

120 

191 

104.105 

130 

192 

104.106 

50 
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193 

105.106 

80 

194 

106.108 

100 

195 

107.109 

90 

196 

101.110 

100 

197 

91  .111 

100 

198 

97  .112 

80  / 

199 

200 


201  EFFORT(I,J,L) 

202  /  1  .20  .AIR  12,  1  .20  .NGS 

10,  1  .20  .SF 

15 

203 

2  .20  .AIR 

14, 

2  .20  .NGS 

8,  2  .20  .SF 

15 

204 

3  .19  .AIR 

14, 

3  .19  .NGS 

10,  3  .19  .SF 

15 

205 

4  .14  .AIR 

6, 

4  .14  .NGS 

12,  4  .14  .SF 

15 

206 

5  .12  .AIR 

4, 

5  .12  .NGS 

9,  5  .12  .SF 

15 

207 

6  .9  .AIR 

2, 

6  .9  .NGS 

15,  6  .9  .SF 

15 

208 

7  .8  .AIR 

12, 

7  .8  .NGS 

6,  7  .8  .SF 

15 

209 

8  .43  .AIR 

10, 

8  .43  .NGS 

6,  8  .43  .SF 

15 

210 

8  .41  .AIR 

14, 

8  .41  .NGS 

18,  8  .41  .SF 

15 

211 

8  .9  .AIR 

12, 

8  .9  .NGS 

15,  8  .9  .SF 

15 

212 

9  .41  .AIR 

16, 

9  .41  .NGS 

15,  9  .41  .SF 

15 

213 

9  .40  .AIR 

18, 

9  .40  .NGS 

15,  9  .40  .SF 

15 

214 

9  .10  .AIR 

12, 

9  .10  .NGS 

9,  9  .10  .SF 

15 

215 

10  .11  .AIR 

16, 

10  .11  .NGS 

10,  10  .11  .SF 

15 

216 

11  .37  .AIR 

14, 

11  .37  .NGS 

12,  11  .37  .SF 

15 

217 

11  .12  .AIR 

16, 

11  .12  .NGS 

28,  11  .12  .SF 

15 

218 

12  .36  .AIR 

4, 

12  .36  .NGS 

29,  12  .36  .SF 

2 

219 

12  .13  .AIR 

6, 

12  .13  .NGS 

24,  12.13.SF 

4 

220 

13  .33  .AIR 

4, 

13  .33  .NGS 

28,  13  .33  .SF 

4 

221 

13  .31  .AIR 

8, 

13  .31  .NGS 

22,  13  .31  .SF 

5 

222 

13  .14  .AIR 

4, 

13  .14  .NGS 

23,  13  .14  .SF 

15 

223 

14  .15  .AIR 

4, 

14  .15  .NGS 

25,  14  .15  .SF 

15 

224 

15  .31  .AIR 

4, 

15  .31  .NGS 

28,  15  .31  .SF 

15 

225 

15  .30  .AIR 

2, 

15  .30  .NGS 

28,  15  .30  .SF 

15 

226 

15  .28  .AIR 

4, 

15  .28  .NGS 

6,  15  .28  .SF 

15 

227 

15  .16  .AIR 

6, 

15  .16  .NGS 

9.  15.16.SF 

2 

228 

16  .23  .AIR 

4, 

16  .23  .NGS 

14,  16  .23  .SF 

5 
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229  16  .17  .AIR 


6,  16  .17  .NGS  15,  16  .17  .SF  4 
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230 

17  .23  .AIR 

231 

17  .22  .AIR 

232 

17  .18  .AIR 

233 

18  .19  .AIR 

234 

19  .22  .AIR 

235 

19  .20  .AIR 

236 

20  .21  .AIR 

237 

21  .22  .AIR 

238 

21  .24  .AIR 

239 

22  .23  .AIR 

240 

23  .28  .AIR 

241 

23  .26  .AIR 

242 

23  .24  .AIR 

243 

24  .25  .AIR 

244 

25  .26  .AIR 

245 

25  .27  .AIR 

246 

25  .68  .AIR 

247 

26  .27  .AIR 

248 

27  .28  .AIR 

249 

27  .62  .AIR 

250 

27  .66  .AIR 

251 

28  .29  .AIR 

252 

29  .61  .AIR 

253 

29  .30  .AIR 

254 

30  .61  .AIR 

255 

30  .60  .AIR 

256 

30  .31  .AIR 

257 

31  .32  .AIR 

258 

32  .33  .AIR 

259 

33  .60  .AIR 

260 

33  .34  .AIR 

261 

34  .38  .AIR 

262 

34  .35  .AIR 

263 

35  .36  .AIR 

264 

35  .37  .AIR 

265 

38  .60  .AIR 

266 

38  .57  .AIR 

267 

38  .39  .AIR 

268 

39  .45  .AIR 

269 

39  .40  .AIR 

270 

41  .42  .AIR 

271 

42  .44  .AIR 

272 

42  .43  .AIR 

2,  17  .23  .NGS 
6,  17  .22  .NGS 
16,  17  .18  .NGS 

10,  18  .19  .NGS 

10,  19  .22  .NGS 
18,  19  .20  .NGS 
14,  20  .21  .NGS 
12,  21  .22  .NGS 
14,  21  .24  .NGS 
12,  22  .23  .NGS 
6,  23  .28  .NGS 
4,  23  .26  .NGS 
6,  23  .24  .NGS 
6,  24  .25  .NGS 
6,  25  .26  .NGS 
6.  25  .27  .NGS 
6,  25  .68  .NGS 
6,  26  .27  .NGS 
6,  27  .28  .NGS 
6,  27  .62  .NGS 
4,  27  .66  .NGS 
4,  28  .29  .NGS 
4,  29  .61  .NGS 
6,  29  .30  .NGS 
4,  30  .61  .NGS 
8,  30  .60  .NGS 
4,  30  .31  .NGS 
4,  31  .32  .NGS 
4,  32  .33  .NGS 
4,  33  .60  .NGS 
4,  33  .34  .NGS 
8,  34  .38  .NGS 
4,  34  .35  .NGS 
2,  35  .36  .NGS 
4,  35  .37  .NGS 
6,  38  .60  .NGS 
4,  38  .57  .NGS 
2,  38  .39  .NGS 
2,  39  .45  .NGS 
2,  39  .40  .NGS 
2,  41  .42  .NGS 
2,  42  .44  .NGS 
2,  42  .43  .NGS 


18. 

17  .23  .SF 

4 

18, 

17  .22  .SF 

4 

18, 

17  .18  .SF 

4 

14, 

18  .19  .SF 

4 

4, 

19  .22  .SF 

2 

8, 

19  .20  .SF 

2 

18, 

20  .21  .SF 

4 

8, 

21  .22  .SF 

2 

8, 

21  .24  .SF 

2 

8, 

22  .23  .SF 

2 

8, 

23  .28  .SF 

2 

8, 

23  .26  .SF 

2 

12, 

23  .24  .SF 

4 

8, 

24  .25  .SF 

2 

8, 

25  .26  .SF 

2 

18, 

25  .27  .SF 

5 

4, 

25  .68  .SF 

2 

4, 

26  .27  .SF 

2 

8, 

27  .28  .SF 

2 

6, 

27  .62  .SF 

2 

10, 

27  .66  .SF 

6 

5, 

28  .29  .SF 

15 

27, 

29  .61  .SF 

15 

25, 

29  .30  .SF 

15 

26, 

30  .61  .SF 

15 

28, 

30  .60  .SF 

15 

22, 

30  .31  .SF 

15 

24, 

31  .32  .SF 

15 

26, 

32  .33  .SF 

15 

24, 

33  .60  .SF 

15 

29, 

33  .34  .SF 

6 

23. 

34  .38  .SF 

2 

28, 

34  .35  .SF 

8 

20, 

35  .36  .SF 

8 

26, 

35  .37  .SF 

2 

27, 

38  .60  .SF 

2 

28, 

38  .57  .SF 

2 

4, 

38  .39  .SF 

2 

9, 

39  .45  .SF 

4 

10, 

39  .40  .SF 

6 

6, 

41  .42  .SF 

2 

10, 

42  .44  .SF 

5 

6, 

42  .43  .SF 

2 

79 


273 

43  .47  .AIR 

2, 

43  .47  .NGS 

274 

44  .45  .AIR 

12, 

44  .45  .NGS 

275 

44  .46  .AIR 

10, 

44  .46  .NGS 

276 

46  .47  .AIR 

16, 

46  .47  .NGS 

277 

47  .48  .AIR 

10, 

47  .48  .NGS 

278 

48  .49  .AIR 

14, 

48  .49  .NGS 

279 

49  .50  .AIR 

16, 

49  .50  .NGS 

280 

49  .51  .AIR 

18, 

49  .51  .NGS 

281 

51  .94  .AIR 

18, 

51  .94  .NGS 

282 

50  .52  .AIR 

18, 

50  .52  .NGS 

283 

52  .94  .AIR 

12, 

52  .94  .NGS 

284 

52  .53  .AIR 

12, 

52  .53  .NGS 
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285 

53  .54  .AIR 

12, 

53  .54  .NGS 

286 

54  .87  .AIR 

2, 

54  .87  .NGS 

287 

54  .83  .AIR 

2, 

54  .83  .NGS 

288 

54  .55  .AIR 

2, 

54  .55  .NGS 

289 

55  .81  .AIR 

4, 

55  .81  .NGS 

290 

55  .82  .AIR 

4, 

55  .82  .NGS 

291 

55  .56  .AIR 

2, 

55  .56  .NGS 

292 

56  .82  .AIR 

4, 

56  .82  .NGS 

293 

56  .73  .AIR 

2, 

56  .73  .NGS 

294 

56  .57  .AIR 

6, 

56  .57  .NGS 

295 

57  .72  .AIR 

6, 

57  .72  .NGS 

296 

57  .59  .AIR 

6, 

57  .59  .NGS 

297 

57  .58  .AIR 

6, 

57  .58  .NGS 

298 

58  .59  .AIR 

4, 

58  .59  .NGS 

299 

59  .61  .AIR 

2, 

59  .61  .NGS 

300 

59  .62  .AIR 

6, 

59  .62  .NGS 

301 

60  .61  .AIR 

4, 

60  .61  .NGS 

302 

62  .63  .AIR 

18, 

62  .63  .NGS 

303 

62  .64  .AIR 

14, 

62  .64  .NGS 

304 

63  .71  .AIR 

14, 

63  .71  .NGS 

305 

63  .64  .AIR 

14, 

63  .64  .NGS 

306 

64  .70  .AIR 

10, 

64  .70  .NGS 

307 

64  .65  .AIR 

10, 

64  .65  .NGS 

308 

65  .70  .AIR 

10, 

65  .70  .NGS 

309 

65  .67  .AIR 

14, 

65  .67  .NGS 

4,  43  .47  .SF  2 

13,  44  .45  .SF  15 

15,  44  .46  .SF  15 

16,  46  .47  .SF  15 

14,  47.48.SF  15 

15,  48  .49  .SF  15 

4,  49.50.SF  15 

20,  49  .51  .SF  15 

8,  51  .94  .SF  2 

4,  50  .52  .SF  2 

10,  52.94.SF  4 

25,  52  .53  .SF  2 
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26,  53  .54  .SF  2 

24,  54  .87  .SF  2 

28,  54  .83  .SF  4 

24,  54  .55  .SF  2 

25,  55  .81  .SF  8 

25,  55  .82  .SF  6 

24,  55  .56  .SF  2 

27,  56.82.SF  2 

20,  56  .73  .SF  4 

22,  56  .57  .SF  5 

26,  57  .72  .SF  2 

24,  57  .59  .SF  2 

26,  57.58.SF  2 

26,  58  .59  .SF  2 

20,  59  .61  .SF  15 

28,  59.62.SF  15 

24,  60  .61  .SF  15 

8,  62  .63  .SF  5 

10,  62  .64  .SF  5 

6,  63  .71  .SF  5 

5,  63  .64  .SF  2 

8,  64  .70  .SF  5 

8,  64  .65  .SF  5 

6,  65  .70  .SF  2 

4,  65  .67  .SF  2 
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310 

65  .66  .AIR 

18, 

65  .66  .NGS 

4, 

65  .66  .SF 

2 

311 

66  .67  .AIR 

14, 

66  .67  .NGS 

10, 

66  .67  .SF 

8 

312 

67  .68  .AIR 

14, 

67  .68  .NGS 

10, 

67  .68  .SF 

15 

313 

67  .69  .AIR 

14, 

67  .69  .NGS 

6, 

67  .69  .SF 

15 

314 

69  .76  .AIR 

18, 

69  .76  .NGS 

8, 

69  .76  .SF 

15 

315 

69  .70  .AIR 

16, 

69  .70  .NGS 

5, 

69  .70  .SF 

15 

316 

70  .71  .AIR 

14, 

70  .71  .NGS 

8, 

70  .71  .SF 

15 

317 

71  .76  .AIR 

18, 

71  .76  .NGS 

8, 

71  .76  .SF 

15 

318 

71  .72  .AIR 

16, 

71  .72  .NGS 

8, 

71  .72  .SF 

15 

319 

72  .75  .AIR 

14, 

72  .75  .NGS 

4, 

72  .75  .SF 

15 

320 

72  .73  .AIR 

18, 

72  .73  .NGS 

18, 

72  .73  .SF 

15 

321 

73  .74  .AIR 

14, 

73  .74  .NGS 

18, 

73  .74  .SF 

15 

322 

74  ,79  .AIR 

16, 

74  .79  .NGS 

18, 

74  .79  .SF 

15 

323 

74  .76  .AIR 

16, 

74  .76  .NGS 

8, 

74  .76  .SF 

15 

324 

74  .75  .AIR 

16, 

74  .75  .NGS 

8, 

74  .75  .SF 

15 

325 

75  .76  .AIR 

14, 

75  .76  .NGS 

4, 

75  .76  .SF 

15 

326 

76  .78  .AIR 

14, 

76  .78  .NGS 

5, 

76  .78  .SF 

15 

327 

76  .77  .AIR 

14, 

76  .77  .NGS 

4, 

76  .77  .SF 

15 

328 

77  .78  .AIR 

12, 

77  .78  .NGS 

8, 

77  .78  .SF 

15 

329 

77  .80  .AIR 

12, 

77  .80  .NGS 

6, 

77  .80  .SF 

15 

330 

78  .79  .AIR 

12, 

78  .79  .NGS 

5, 

78  .79  .SF 

15 

331 

79  .81  .AIR 

14, 

79  .81  .NGS 

18, 

79  .81  .SF 

15 

332 

79  .80  .AIR 

12, 

79  .80  .NGS 

10, 

79  .80  .SF 

5 

333 

80  .81  .AIR 

14, 

80  .81  .NGS 

4, 

80  .81  .SF 

2 

334 

80  .83  .AIR 

10, 

80  .83  .NGS 

4, 

80  .83  .SF 

3 

335 

80  .84  .AIR 

12, 

80  .84  .NGS 

8, 

80  .84  .SF 

2 

336 

81  .83  .AIR 

2, 

81  .83  .NGS 

8, 

81  .83  .SF 

2 

337 

81  .82  .AIR 

2, 

81  .82  .NGS 

8, 

81  .82  .SF 

5 

338 

83  .86  .AIR 

2, 

83  .86  .NGS 

4, 

83  .86  .SF 

2 

339 

83  .84  .AIR 

2, 

83  .84  .NGS 

10, 

83  .84  .SF 

6 

340 

84  .85  .AIR 

2, 

84  .85  .NGS 

8, 

84  .85  .SF 

2 
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341 

85  .106. AIR 

4. 

85  .106.NGS 

8, 

85  .106.SF 

15 

342 

85  .86  .AIR 

4, 

85  .86  .NGS 

8, 

85  .86  .SF 

15 

343 

86  .89  .AIR 

8, 

86  .89  .NGS 

4. 

86  .89  .SF 

15 

344 

86  .87  .AIR 

8, 

86  .87  .NGS 

20, 

86  .87  .SF 

15 

345 

87  .88  .AIR 

6, 

87  .88  .NGS 

24. 

87  .88  .SF 

15 

346 

88  .89  .AIR 

6, 

88  .89  .NGS 

4, 

88  .89  .SF 

15 

347 

89  .103. AIR 

4, 

89  .103.NGS 

4, 

89  .103.SF 

15 
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348 

89  .100. AIR 

2, 

89  .100.NGS 

13,  89.100.SF 

15 

349 

89  .90  .AIR 

6, 

89  .90  .NGS 

4,  89  .90  .SF 

15 

350 

90  .91  .AIR 

4, 

90  .91  .NGS 

5,  90  .91  .SF 

15 

351 

91  .99  .AIR 

18, 

91  .99  .NGS 

11,  91  .99  .SF 

15 

352 

91  .98  .AIR 

16, 

91  .98  .NGS 

5,  91  .98  .SF 

15 

353 

91  .97  .AIR 

14, 

91  .97  .NGS 

4,  91  .97  .SF 

15 

354 

91  .96  .AIR 

14, 

91  .96  .NGS 

12,  91  .96  .SF 

15 

355 

91  .92  .AIR 

10, 

91  .92  .NGS 

12,  91  .92  .SF 

15 

356 

92  .93  .AIR 

10, 

92  .93  .NGS 

12,  92.93.SF 

15 

357 

93  .96  .AIR 

14, 

93  .96  .NGS 

4,  93  .96  .SF 

15 

358 

93  .94  .AIR 

12, 

93  .94  .NGS 

4,  93.94.SF 

15 

359 

94  .95  .AIR 

12, 

94  .95  .NGS 

8,  94  .95  .SF 

15 

360 

95  .97  .AIR 

12, 

95  .97  .NGS 

10,  95  .97  .SF 

15 

361 

95  .96  .AIR 

12, 

95  .96  .NGS 

10,  95  .96  .SF 

15 

362 

98  .99  .AIR 

14, 

98  .99  .NGS 

8,  98.99.SF 

15 

363 

99  .lOl.AIR 

12, 

99  .101. NGS 

10.  99  .lOl.SF 

15 

364 

99  .102. AIR 

12, 

99  .102.NGS 

8,  99  .102.SF 

15 

365 

99  .100. AIR 

12, 

99  .lOO.NGS 

6,  99  .lOO.SF 

15 

366 

102. 103.  AIR 

12, 

102.103.NGS 

6,  102.103.SF 

15 

367 

103. 107.  AIR 

12, 

103.107.NGS 

10,  103.107.SF 

15 

368 

103.104.AIR 

10, 

103.104.NGS 

14,  103.104.SF 

15 

369 

104.107.AIR 

10, 

104.107.NGS 

12,  104.107.SF 

15 

370 

104.105.AIR 

16, 

104.105.NGS 

12,  104.105.SF 

16 

371 

104.106.AIR 

10, 

104.106.NGS 

8,  104.106.SF 

14 

372 

105.106.AIR 

10, 

105.106.NGS 

8,  105.106.SF 

14 

373 

106. 108.  AIR 

14, 

106.108.NGS 

10,  106.108.SF 

15 

374 

107.109.AIR 

14, 

107.109.NGS 

8,  107.109.SF 

14 

375 

lOl.llO.AIR 

12, 

lOl.llO.NGS 

10,  101.1  lO.SF 

16 

376 

91  .111.AIR 

12, 

91  .lll.NGS 

10,  91  .lll.SF 

16 

377 

97  ii2.AIR 

12, 

97  .112.NGS 

8,  97  .112.SF 

14/ 

378 

379  ASSET  (L) 

380  /  AIR  4 

381  NGS  4 

382  SF  4  /  ; 

383 

384  POSITIVE  VARIABLE 

385  H(I,J); 

386 

387  VARIABLE 

388  A(I) 

389  MAXCAP; 

390 
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391  BINARY  VARIABLES 


392 

G(I,J,L); 

393 

394 

A.FX(T’) 

=  0; 

395 

A.FX(’2’) 

=  0; 

396 

A.FX(’3’) 

=  0; 
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397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 


A.FX(’4’)  =  0; 
A.FXC5’)  =  0; 
A.FX(’6’)  =  0; 
A.FX(’7-)  =  0; 
A.FX(’108’)  =  1; 
A.FX(’109’)  =  1; 
A.FX(’110’)  =  1; 
A.FX(’lir)  =  1; 
A.FX(’112’)  =  1; 


EQUATIONS  OBJ 
ARCl(U) 
ARC2(I,J) 
WPN(L) 


define  objective  function 
equation  for  forward  arc 
equation  for  backward  arc 
weapon  expenditure  for  each  arc; 


•  >  >  >  minimize  <  <  < 

OBJ.. 

MAXCAP  =E=  SUM  (  (U)  $  (  CAP(I,J)  GT  0  )  ,  CAP(U)  • 
H(U)); 


*  >  >  >  subject  to  <  <  < 


ARC1(I,J)  $  (  CAP(I,J)  GT  0).. 

A(I)  -  A(J)  +  H(I,J)  +  SUM  (L,  G(U,L))  =G=  0; 


ARC2(I,J)  $  (  CAP(I,J)  GT  0).. 

A(J)  -  A(I)  +  H(I,J)  +  SUM  (L,  G(I,J,L))  =G=  0; 

WPN(L).. 


83 


428  SUM(  (I,J)  $  (  CAP(I,J)  GT  0  ),EFFORT(I,J,L)  *  G(I,J,L))  =L= 

429  ASSET(L)  ; 

430 

431  MODEL  NETINT  /ALL/; 

432  SOLVE  NETINT  USING  MIP  MINIMIZING  MAXCAP; 

433  DISPLAY  MAXCAP.L; 

434  DISPLAY  H.L; 

435  OPTION  G:0;2;l;  DISPLAY  G.L; 


COMPILATION  TIME  =  0.380  SECONDS 
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MODEL  STATISTICS  SOLVE  NETINT  USING  MIP  FROM  LINE  432 

GAMS  2.19  IBM  CMS 


MODEL  STATISTICS 


BLOCKS  OF  EQUATIONS  4 
BLOCKS  OF  VARIABLES  4 
NON  ZERO  ELEMENTS  2817 


SINGLE  EQUATIONS  356 
SINGLE  VARIABLES  817 
DISCRETE  VARIABLES  528 


GENERATION  TIME  =  1.550  SECONDS 


EXECUTION  TIME  =  1.630  SECONDS 

1 
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SOLUTION  REPORT  SOLVE  NETINT  USING  MIP  FROM  LINE  432 

GAMS  2.19  IBM  CMS 


SOLVE  SUMAR Y 


MODEL  NETINT 
TYPE  MIP 
SOLVER  ZOOM 


OBJECTIVE  MAXCAP 
DIRECTION  MINIMIZE 
FROM  LINE  432 


\ 


••••  SOLVER  STATUS  1  NORMAL  COMPLETION 
•••*  MODEL  STATUS  1  OPTIMAL 
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•  •  ^  •  OBJECTIVE  VALUE  20.0000 

RESOURCE  USAGE,  LIMIT  1.162  10000.000 

ITERATION  COUNT,  LIMIT  303  90000 

Z  O  M  /  X  M  P  ---  Version  2.1  Jun  1988 

Courtesy  of  Dr  Roy  E.  Marsten, 

Department  of  Management  Information  Systems, 
University  of  Arizona, 

Tucson  Arizona  85721,  U.S.A. 


PROBLEM  SPECIFICATIONS 


BEGIN 


•  SPECS  FILE,  VERSION  2.1  JUN  1987 

« 

•  AMOUNT  OF  PRINTOUT  IN  STATUS  FILE 

PRINT  CONTINUOUS  0 
PRINT  BRANCH  0 
PRINT  HEURISTIC  0 
PRINT  TOUR  0 


PARAMETERS  CONTROLLING  LP 


•  PARAMETERS  CONTROLLING  THE  HEURISTIC 

« 

HEURISTIC  YES 

•HEURISTIC  NO 

« 

•  PARAMETERS  CONTROLLING  THE  BRANCH  AND  BOUND. 

• 

BRANCH  YES 
QUIT  NO 
DIVE  YES 
INCUMBANT= 300.0 

INVERT  =  50 

« 

MAX  SAVE  5 
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EXPAND  3 

SELECT  3 

» 

08/04/91 
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SOLUTION  REPORT  SOLVE  NETINT  USING  MIP  FROM  LINE  432 

GAMS  2.19  IBM  CMS 


END 

Work  space  needed  (estimate)  --  55974  words. 
Work  space  available  --  55974  words. 

Maximum  obtainable  -  295165  words. 

The  LU  factors  occupied  1319  slots  (estimate  7388). 

Iterations:  Initial  LP  303,  Time:  0.83 

Heuristic  0,  0.00 

Branch  and  bound  0,  0,00 

Final  LP  0,  0.00 


••••  REPORT  SUMMARY  :  0  NONOPT 

0  INFEASIBLE 
0  UNBOUNDED 

08/04/91 
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EXECUTING 
GAMS  2.19  IBM  CMS 


--  433  VARIABLE  MAXCAP.L  =  20.00 


— -  434  VARIABLE  H.L 

52 

50  1.00 


435  VARIABLE  G.L 
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AIR  NGS  SF 


51  .94 

53  .54 

54  .87 

83  .86 

84  .85 


1 

1 


•••*  FILE  SUMMARY  FOR  USER  8847P 

INPUT  SC2  GAMS  A 
OUTPUT  SC2  LISTING  A 

EXECUTION  TIME  =  0.270  SECONDS 
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